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1.  INTRODUCTION 

The  ABORC  (arbitrary  body  of  revolution)  computer  code  for  the  solu- 
tion of  system-generated  electromagnetic  pulse  (SGEMP)  and  internal  elec- 
tromagnetic pulse  (lEMP)  on  arbitrary  bodies  of  revolution  is  documented 
in  this  report.  The  code  is  designed  primarily  for  tlie  solution  of  elec- 
tromagnetic currents  and  fields  produced  by  arbitrary  axisymmetric  emis- 
sion of  electrons  due  to  photons  incident  on  objects  such  as  satellites 
or  missiles.  The  complete  set  of  Maxwell's  equations  is  solved  and 
coupled  with  particle  motion  representing  currents.  Self-consistent, 
time-dependent  solutions  of  currents  and  fields  are  obtained  in  two 
dimensions . 

This  report  contains  a description  of  code  capabilities  and  various 
tests  which  have  been  performed  to  determine  the  validity  of  its  solu- 
tions by  comparing  ABORC  fields  and  currents  with  analytical  solutions, 
where  available,  and  also  with  other  computer  codes. 

Detailed  descriptions  of  the  physics  and  modeling  are  contained  in 
Appendix  A.  Sample  ABORC  calculations  on  problems  of  interest  in  SGEMP 
are  found  in  Appendices  B and  C.  Results  in  Appendix  B permit  ordor-of- 
magnitude  estimates  for  SGEMP  responses  over  wide  ranges  of  object  size 
and  photon  excitation.  Af)pcndix  C contains  pertinent  considerations  of 
geometry  effects  under  space-charge- 1 imited  (SCL)  conditions.  The  appen- 
dices show  by  illustration  the  range  of  applicability  of  ABORC  to  SCiEMP 
ca 1 cul at i ons . 

Some  numerical  sensitivity  considerations  are  discussed  in  Appendix 
D.  These  are  designed  to  aid  the  user  in  ju'oducing  results  containing 
minimal  statistical  noise.  A complete  user's  manual  for  ABORC  and  related 
gra()hics  is  jirovided  in  Appei.dix  li. 


2.  ABORC  COMPUTER  CODE  DESCRIPTION 


The  ABORC  code  solves  Maxwell's  equations  with  self-consistent  elec 
tron  motion  in  a vacuum  for  axisymmetric  geometries.  Direct  finite- 
differencing  of  the  field  equations  is  done  employing  generalized  coordi 
nates,  and  finite  "particles"  charge  are  followed  through  the  spatial 
mesh  of  zones  to  obtain  currents.  Emission  of  arbitrary  energy,  angular 
spatial,  and  time  distributions  of  currents  can  be  specified,  including 
fully  time-dependent  spectrum  and  spatial  distributions.  Randomizing 
techniques  are  employed  for  all  distributions  for  efficient  numerical 
representation. 

Boundary  conditions  currently  available  in  the  code  require  the 
specification  of  an  outer,  perfectly  conducting  cylinder.  Free-space 
solutions  can  be  obtained  by  moving  the  outer  boundary  out  so  the  clear 
time  (the  time  in  which  reflections  from  the  outer  wall  return  to  the 
structure]  is  larger  than  the  problem  time  of  interest.  The  shape  of 
the  inner  wall  of  the  outer  conductor  can  be  modified  to  an  arbitrary 
body  of  revolution  by  specifying  conductivities  within  the  outermost 
cylinder,  but  the  final  spatial  extent  of  the  calculation  is  the  above- 
mentioned  perfectly  conducting  cylinder. 

Finite  conductivities  can  be  specified  representing  imperfect  con- 
ductors, and  dielectric  structures  may  be  treated  by  specifying  proper 
dielectric  constants.  While  the  former  is  programmed  into  ABORC  and  is 
implem.ented  by  simple  input  card  specification,  the  latter  capability 
is  somewhat  limited.  Dielectrics  with  e = (free  sjiace  permittivity) 
can  be  s{)ecified  throughout  the  volume  of  the  calculation  by  simply 
specifying  conductors  of  very  low  conductors.  Resistors  are  specified 
in  this  way,  where  t!ie  conductivity  is  chosen  to  give  the  sought-after 
resistance  value.  Care  must  also  be  taken  liere  to  ensure  that  the  skin 
depth  of  the  resistor  is  large  enough  to  allow  proper  field  penetration. 
Dielectrics  witti  c /-  e.^  pi'esent  a different  problem,  however,  in  that  no 


provision  has  been  made  for  them  in  the  field  equations.  Programming 
changes  are  therefore  necessary  to  add  this  capability.  For  certain 
ca.ses  with  small,  well  defined  volumes  of  dielectrics,  'nowever,  this 
change  can  be  relatively  minor. 

While  most  SGF.MP/IBMP  calculations  are  performed  assuming  complete 
absorption  of  photo-electrons  upon  impact  with  surfaces,  back-scattering 
of  electrons  can  be  specified  in  ABORC  if  desired.  Variable  amounts  of 
charge  with  modified  energy  and  direction  of  propagation  can  be  re-emitted 
from  surfaces  upon  contact.  This  feature  is  presently  limited  to  back- 
scattering  from  a simple  cylinder  with  energy-  and  angle-independent 
reflection  fractions,  energy  dissipation,  and  specular  angle.  Studies 
undertaken  with  this  code  option  have  shown  considerable  effect  on  SGEMP 
response  under  certain  conditions  (Ref.  1) . 

Graphics  features  include  the  capability  of  simple  specification  of 
most  major  calculational  quantities  at  arbitrary  locations  to  be  dis- 
posed to  various  plotting  routines.  These  routines  include  printer 
plots,  pen  plots,  CRT  plots,  and  16mm  computer-generated  movies  of  elec- 
tron trajectories.  Files  of  data  from  different  ABORC  runs  can  be  over- 
laid conveniently,  allowing  accurate  analysis  of  calculations.  Time 
histories  and  spatial  distributions  of  fields  and  currents  can  be  corre- 
lated directly  with  particle  trajectories  using  the  movie  capability. 

The  movies  provide  a helpful  tool  in  debugging  complicated  geometry  cal- 
culations where  electron;'  arc  emitted  from  many  surfaces  and  spaco-cliargc- 
li mi  ting  may  occur. 


, 1’.  Weiiaas  and  A.  . Wootis,  ''Fdeex  roii  llackst.if  x e r i ng  iltt'cts  on 
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I lad.  Con  fe  rtiice  on  Nnclx'.ir  f,  Space  Radi. it  ion  I'Cffi^ts,  du')' 


5 


3.  ABORC  VERIFICATION 


ABORC  has  been  checked  against  known  solutions  of  Maxwell's  equations 
and  against  other  computer  codes.  The  results  of  the  checkouts  are  dis- 
cussed here.  Comparisons  are  for  both  non-space-charge-limited  and  space- 
charge-limited  currents. 


3.1  EMPTY  CYLINDER  FIELD  SOLUTION 


ABORC  has  been  compared  with  analytic  solutions  for  electric  fields 
in  an  empty  cylinder.  These  analytic  solutions  can  be  derived  when  a sinu- 
soidal axial  current  density  having  a Bessel  function  radial  dependence  is 
specified: 


for  0 < u < L,  zero  otherwise,  (1) 


where 

R = cylinder  radius, 

L = cylinder  length, 
r,z  = radial,  axial  coordinates, 
u = vt  - z , 

V = phase  velocity, 

= first  zero  of  Bessel  function  .1^^. 

i'hc  solutions  for  the  radial  and  axial  electric  fields  at  time  = L/V 


0 


(3) 


A test  problem  was  run  with  ABORC  where  the  axial  current  was  speci- 
fied according  to  Eq.  1.  The  test  problem  had  the  following  characteristics. 
Cylinder  radius  (R)  = 0.5  m 
Cylinder  length  (L)  = 1.0  m 
Phase  velocity  (v)  = 0.2c 

These  parameter  inputs  resulted  in  an  electron  flight  time  of  16.7  nsec 
across  the  cavity,  which  was  used  as  the  maximum  running  time. 

ABORC  was  run  with  both  constant  and  variable  spatial  zones  as  a 
check  of  the  variable-zone  capability  of  the  code.  The  following  inputs 
pertaining  to  grid  sizes  were  employed. 

Number  of  radial  zones:  10 
Number  of  axial  zones:  20 

Time  step:  2.09  x 10  sec 

Zone  sizes  varied  by  up  to  a factor  of  10  in  the  variable-zoned  case. 

Equation  1 was  used  to  update  the  axial  current  density  each  time 
step.  Fields  from  the  code  were  compared  with  the  analytic  solutions  (Eqs. 

2 and  3)  at  16.7  nsec.  The  results  of  the  comparison.s  are  seen  in  Figures 
1 and  2 for  the  axial  electric  field  at  the  emitting  end  of  the  cylinder 
and  the  radial  electric  field  at  the  side  of  the  cylinder,  respectively. 

I’he  arrows  on  the  geometry  figures  indicate  tlie  electron  direction  of 
propagat ion . 

AB0R(!  results  are  from  the  varialile-zone  case,  which  agieeii  well 
with  the  constant-zone  case.  Agreement  of  tlie  code  and  analytic  results 
is  very  gooii.  SI  iglit  differences  may  tie  attributed  to  numerical  inaccura- 
cies resulting  from  finite  grid  spacings.  I'his  test  of  the  code  iiulicatcs 
that  the  field  calculation  portion  is  operating  correctly,  given  the  cur- 
in  an  empty  cyl inder. 


rent 


son  of  \RORr  with  analytic  solution  for  empty  cylinder;  radial 
c field  at  cylinder  wall  versus  axial  position  at  time  = 16.67 
r va ri able- increment  spatial  toning 


3.2  DOUBLE-CYLINDER  FIELD  SOLUTION 


ABORC  has  been  compared  with  analytic  solutions  for  E-  and  H-fields 
in  a double-cylinder  geometry  (Figure  3).  Analytical  solutions  to  Max- 
well's equations  are  particularly  easy  to  determine  for  the  symmetry  and 
initial  conditions  considered  here.  This  simplification  follows  from  the 
fact  that  the  H-field  in  the  azimuthal  direction  satisfies  the  boundary 
condition  on  a perfect  conductor  automatically.  Any  specification  of  the 
E-field  which  satisfies  the  boundary  conditions  and  is  zero  initially 
implies  a value  of  the  magnetic  field  from  the  curl  E equation.  Tiie  com- 
bination of  E-  and  H-fields  then  determines  the  value  of  the  current  den- 
sity from  the  curl  H equation,  and  it  is  that  current  which  must  be  used 
to  drive  the  computer  solution  to  produce  the  E-  and  H-fields. 
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An  example  of  these  solutions  is  shown  in  the  following  test  problem. 
Analytic  solutions  are  known  by  the  above  methods  when  an  axial  current 
density  having  the  following  dependence  is  specified. 

■^z  " ^'^i  ^0  ^*^i^  ~ *^0  ^^i  (t) 

for  0 < u < L,  zero  otherwise, 

where 


k.  = 3.1228  for  b/a  - 2, 

1 

b = outer  cylinder  radius, 
a = inner  cylinder  radius, 

L = cylinder  length, 
r,z  = radial,  axial  coordinates, 
u = vt  - z , 

V = phase  velocity, 

are  Bessel  functions  of  order  zero. 

0 0 


The  value  of  k^  causes  the  quantity  in  brackets  to  vanish  at  r = a and 
r = b.  The  axial  and  radial  electric  fields  at  a time  equal  to  the  flight 
time  for  electrons  down  the  length  of  the  cylinders  are  given  by 


E 

z 


-TT  1 
L-  e^v 


Ikp  jQ(k.  r/al  - Yp(k.  r/a) 


cos 


TTZ 

L 


(5) 


and 
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A test  problem  was  run  with  ABOUC  where  the  axial  current  density  was 
specified  according  to  bq.  4.  'Die  problem  had  the  following  characteristics. 


Inner  cylinder  radius,  a. 

0.25  m 

Outer  cylindt-r  radius,  b: 

O.S  11! 

(iylindi.r  length,  1.: 

1 m 

i’eak  cu.rrenl'  density, 

'cO , Kfl 

(occurs  at  a time  of  one-- 

b.alf  an 

elc'.tron  flight:  time  at  a rrrdial 
piisit  ion  midway  between  c\linders' 
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Phase  velocity,  v: 

0.2c 

Number  of  radial  zones  between  cylinders: 

10 

Number  of  axial  zones: 

20 

Constant  spatial  zone  sizes: 

Ar  - 

0.25 

m 

Az  = 

0.05 

m 

Time  step: 

4 . 17 

.X  10' 

-11 

The  problem  characteristics  resulted  in  an  electron  flight  time  of  16.7  nsec 
The  electric  fields  normal  to  the  inner  and  outer  cylinder  walls  and 
the  end  away  from  the  emission  surface  are  plotted  at  the  time  of  16.7  nsec 
in  Figures  4,  5,  and  6.  The  arrows  on  the  geometry  figures  indicate  the 
election  direction  of  propagation.  The  analytic  curves  arc  also  shown. 
Agreement  of  the  code  and  analytical  results  is  excellent,  with  minor  dif- 
fex  aces  attributable  to  finite  grid  sizes  used  in  ABORC.  This  test  indi- 
cates that  the  field  solution  is  operating  coi,.ectly  in  cylinders  with 
objects  inside. 


3.3  PARTlChL;  P'MISSION  TliSTR 

Several  checkouts  of  the  particle-emission  portion  of  the  \BORC  code 
have  been  made.  Some  of  the  tests  are  described  here.  The  particle  emit- 
ter has  built-in  coding  to  store  information  on  each  emission  energy  dis- 
tribution used  in  a given  calculation.  This  information  can  be  plotted  at 
the  end  of  a run  and  compared  with  the  desired  input  energ)  distributions 
to  test  statistics.  ro’’‘pari sons  of  these  "desired"  and  "obtained"  distri- 
butions have  shown  agreements  witliin  iJeniiissible  statistical  deviations. 

The  deviations  are  due  to  the  Monte  Carlo  techniques  cmplcycd  by  the  code. 

Emissions  of  currents  from  a variety  of  geometric  objects  have  been 
studied.  In  many  cases,  the  total  charge  emitted  was  analyti.:al  ly  calcu- 
lable (such  as  emission  from  simple  cylinders  or  cones).  The  emitter  has 
been  verified  under  these  cempari sons , with  pulses  botii  long  and  short 
compared  to  object  dimensions.  The  latter  case  provides  a test  of  the 
delayed  emission  capability  of  the  co  e. 

A particularly  complete  eh;  ckout  of  the  particle  emitter  was  afforded 
by  modeling  a time-dependent  spectrum  in  an  elect  ron-t)eam  simulation  cal- 
culation, Tile  beam  was  modeled  with  a triangular  pulse  shajie,  and  the 


energy  of  the  electrons  was  specified  as  triangular  in  time.  These  sim- 
plified representations  permitted  an  analytical  solution  to  he  obtained 
for  the  time-averaged  beam  energy.  This  energy  agrees  with  that  obtained 
from  ABORC  when  the  total  electron  kinetic  energy  emitted  is  divided  by 
the  total  charge  emitted.  The  latter  numbers  are  standard  ABORC  printout 
quantities . 

3.4  CHECKOUTS  INCLUDING  PARTICLE  MOTION 

Several  checkouts  of  the  particle  motion  of  ABORC  have  been  made 
against  analytic  solutions  and  existing  computer  codes  in  both  non-space- 
charge-  1 imited  and  space-charge-limited  situations. 

3.4.1  Anal ytical  Compari sons  for  Par tic le  Motion 

Behavior  of  individual  particles  injected  into  an  empty  cylinder  has 
been  examined  under  field  conditions  where  analytic  solutions  are  deriv- 
able. In  three  difierent  tests,  each  of  the  three  field  quantities  — 

E^,  and  --  have  been  specified  constant  throughout  the  cavfty,  with  the 
other  two  being  set  to  zero.  Then  particles  were  injected  in  such  a way 
that  simple  analytical  expressions  could  be  derived  to  compare  with 
resultant  particle  motion.  The  particles  were  injected  parallel  to  E_^ 
and  and  along  the  axis  when  was  non-zero.  The  particles  wore  found 
to  come  to  a stop  at  positions  within  5%  of  the  analytic  predictions  when 
about  20  time  steps  from  the  initial  to  final  positions  were  taken.  The 
radius  of  curvature  for  particles  injected  perpendicular  to  the  l!-field 
was  found  to  agree  very  well  with  the  analytical  value,  slight  differences 
being  attributed  to  finite  gria  sizes. 

Obviously,  the  abovementioned  particle  checkouts  are  limited  in  scope, 
but  they  do  provide  some  confidence  in  results.  More  comprehensive  tests 
are  rejiorted  below. 
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3.4.2  Comparisons  of  ABORC  with  Other  SGliMP  Codes 

ABORC  has  been  compared  with  the  computer  codes  DYN'ACYI.  (l\efs.  2-4) 
and  SliMP  (Refs.  .S,6)  for  the  configuration  of  a cylinder  within  a cylinder. 
Both  of  these  codes  solve  Maxwell's  equations  employing  self-consistent 
particle  motion  for  describing  photo-electron  currents.  The  codes  have 
been  applied  to  varieties  of  ICMP  and  SGEMP  problems.  These  comparisons 
compare  the  ability  of  ABORC  to  treat  an  object  interior  to  the  enclo- 
sure and  also  to  treat  space-charge- 1 imi ted  situations. 


Non-Space-Charge- Limited  Case  - DYNACYL  Compari son 


The  test  problem  has  the  following  characteristic  dimensions. 


Geometry 

Outer  cylinder  length; 
Outer  cylinder  diameter: 
Inner  cylinder  length: 
Inner  cylinder  diam.eter: 
Inner  cylinder  position: 
Emission  current 
Pulse  shape: 

Peak  current; 

Spatial  distribution: 

Electron  velocity 
Angular  distribution; 
Space-charge- 1 i mi  t ing  did  not  occi.r 


6.0  m 
6.0  m 
0.8  m 
1 . 4 m 

center  of  outer  cylinder 

triangular  (symmetric),  20-nscc  FlVliM 
1 amp 

electrons  emitted  from  top  of  inner 
cyl inder  only 

0.2c,  axial  direction 

stra ight-out  cm i ss i on 

here  because  the  fields  were  not 


allowed  to  affect  electron  motion  in  this  test. 


‘-'1.  N.  belnier  ct  al.,  "SGEMI’’  Plienotnenology  and  Computer  Code  Devel- 
opment," DNl  36.S3F,  November  11,  l'.)74. 

‘"’E.  P.  Wenaas  and  A.  .1.  Woods,  "CoiniKir  i sons  of  Quasi -Static  and  Fully 
Dvnamic  Solutions  for  Electromagnetic  Field  Calculations  in  a Cylindrical 
(iavity,"  1 EFE  T rans . Nuc  1 . Sc  i . NSj^2J^,  December  1074,  j).  2.30. 

'^D.  C.  Osbcni  et  a!.,  "Large-Area  L lect  ron- Bearn  L.xpei'iment  s , " INi'FL- 
RT  8 11)1-011,  .July  l.S,  107.3. 

■’D.  L.  Mang;in  ;ind  K.  A.  Perala,  "Satellite  SCLMl’  Suri'ace  Current  i'ecii - 
nicjiics,"  IFLL  Frans.  Nucl.  Sci.  N.S-22,  No.  (> , December  107.3,  [i.  2420. 

. Stettner  and  II.  L.  Longley,  "Description  of  the  SCliMP  ComjMiter 
Code  SLMP,"  Mission  Research  Corp.  rejHrrt  MRC-R-ldl,  .June  1073. 


0.25  nsec,  ABORC 
0.33  nsec,  DYNACYL 


Grid  sizes  employed  in  the  two  codes  are  listed  below.  Constant 
spatial  zoning,  equal  in  each  direction,  was  used  in  both  cases. 

Number  of  radial  zones:  15 

Number  of  axial  zones:  30 

Maxwell's  equation 
time  step 

The  magnetic  field  near  the  edge  of  the  inner  cylinder  farthest  from 
the  emission  point  is  shown  for  the  two  codes  in  Figure  7.  This  field 
position  is  on  the  inner  cylinder  v/all  at  a radial  position  of  0.7  m.  The 
axial  position  is  about  0.2  m from  the  end  of  the  inner  cylinder.  The 
curves  show  good  overall  agreement,  with  differences  in  the  start-up  times 
attributable  to  emission  current  pulse  start-up  time  differences  between 
the  two  cases.  These  differences  were  due  to  initial  electron  position 
and  time  step  treatments  existing  in  the  two  codes  at  the  time  of  this 
test,  and  also  to  field  point  position  differences  due  to  finite  zone 
sizes. 


Space-Charge-Limited  Case  - SEMP  Comparison 

ARORC  has  been  applied  to  a space-charge-limited  calculation  of  fields 
and  currents  in  a double-cylinder  geometry  and  results  compared  with  those 
obtained  by  Higgins  for  similar  problem  conditions.  Unfortunately  for 
this  comparison,  the  emission  current  spectra  differed  slightly,  but  all 
other  problem  characteristics  were  approximately  the  same.  The  input 
conditions,  while  not  perfect,  do  permit  comparison  of  results  of  testing 
most  facets  of  the  codes,  including  the  ability  to  perform  space-charge- 
limited  calculations. 

Problem  conditions  are  listed  in  Table  1.  These  and  the  results  are 

cast  in  the  scaling  law  forms  discussed  in  Ai^pendix  B.  These  jjroblem  con- 

- 2 

ditions  result  in  a peak  surface  current  of  about  1.6  x 10  " times  the 
peak  emission  current,  indicating  a highly  sjiacc-chargc- 1 imi ted  situation. 

Results  for  surface  currents  at  two  positions  on  the  cylinder  are 
given  in  Figures  8 and  9 for  the  two  codes.  The  curves  slww  fairly  good 
agreement,  with  results  becoming  more  similar  as  the  field  poin.t  is  far- 
ther removed  from  the  [)article  region.  SliMP  currents  appear  to  lie  some- 
what more  space-charge- 1 imi ted  than  ABORC  results.  This  is  consistent 

*1).  Higgins,  MRC  Corp.,  [irivate  communication,  March  197(). 
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Table  1 

ABORC/SEMP  COMPARISON  PROBLEM 


Geometry: 


Emission  Current; 
Pulse  shape; 
Pulse  rise  time; 


Level  ; 

Spatial  distribution: 


Double  cylinders,  length  = diameter 

Outer  diameter  = 3X  inner  diameter 

Inner  cylinder  connected  by  wire  to 
outer  cylinder 


Triangular,  symmetric 

T = 0.25  X (2Tir/c),  where 

r = inner  cylinder  radius, 

T = pulse  rise  time 
c = 3 X 10®  m/sec 

5 . S X 10^  amp 

Uniform  emission  over  end  of  inner 
cylinder  away  from  wire 


Angular  distribution: 


Cos  0 measured  from  surface  normal 


.8 


.5 


2. 


ABOiU'  ail' 
r a cyl  iiidei 


SliMl’  surface 
for  '^;iace-i  iiar 


with  observed  emission  electron  spectrum  differences  where  the  distri- 
bution used  in  SEMP  had  slightly  lower  average  energy  than  that  used  in 
ABORC.  Response  of  the  object  to  the  reflected  wave  from  the  outer  walls 
is  seen  clearly  in  Figure  8 in  both  code  results.  The  clear  time  is 
indicated  on  the  graphs.  In  general,  the  codes  show  enough  agreement  to 
suggest  that  both  are  giving  reasonable  estimates  for  space-charge- 
limited  surface  currents.  Further  investigation  of  the  differences 
between  the  results  would  require  analyzing  the  details  of  the  problem 
numerics  and  code  formulations. 

3.4.3  Summary  of  Code  Verifications 

A summary  of  the  code  checkouts  is  given  below. 

Empty-cylinder  field  solution  - analytic 
Double-cylinder  field  solution  - analytic 
Particle-emitter  tests  - analytic 
Particle  motion  tests  - analytic 
Comparison  witli  DYNACYi,  - non-SCL 
Comparison  with  SEMP  - SCL 

In  addition  to  these  checkouts,  ABORC  results  were  included  in  a compar- 
ison of  data  From  several  com])uter  codes  and  analytical  predictions  con- 
ducted by  del’lomb  (Ref.  7).  Problem  conditions  varied  slightly  between 
ABORC  calculations  and  tlie  other  models,  but  cominirison  of  ihe  data  was 
still  a meaningful  exercise.  ABORC  surface  current  I'leji)'  values  and  times 
of  occurrence  were  found  to  he  in  g.ood  general  agreement  with  the'.e  other 
sources  over  wide  ranges  of  emission  currents  and  specta  including  space- 
c ha  rg e - 1 i m i t ed  c ond i t i on  s . 


7 

' 1: . P.  ilel’lomb,  "Analytical  Predictions  of  Sfil.Ml’  Response  and  tioni 
pai'isons  with  tiompnter  (f-i  1 cu  1 at  i ons  , " K.aman  Sciences  iwpurt 
to  he  published  A|ii-il  197(). 
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4 . SUMf4ARY 

The  ABORC  computer  code  for  the  solution  of  SGEMP/IEMP  problems  for 
arbitrary  bodies  of  revolution  in  two  dimensions  has  been  documented  in 
this  report.  A description  of  code  capabilities  and  code  verifications 
are  given  in  the  main  body,  while  contents  of  appendices  include  a 
detailed  description  of  ABORC  physics  and  modeling,  sample  calculations, 
a user's  manual,  and  numerical  grid  choice  considerations. 

A brief  description  summarizing  ABORC  capabilities,  physics,  and 
modeling  is  given  in  Table  2.  ABORC  verification  included  testing  vir- 
tually every  code  capability  to  some  degree,  with  at  least  reasonable 
similarity  to  other  solutions  being  obtained.  Tested  code  portions 
include  empty-cylinder  fields,  non-empty-cylinder  fields,  individual 
particle  trajections,  and  botli  non-.SCL  and  SCI.  SCIiMP  calculations  in 
double-cylinder  geometries.  ABORC  has  also  been  compared  with  DYNASPIIERE 
(Ref.  2]  on  a .spherical  geometry  under  SCI.  conditions,  with  reasonable 
results.  Thus,  the  code  appears  to  have  the  ca])ability  to  handle  curved 
surfaces  as  well.  Verification  of  ABORC  on  very  com]')l  icated  geometries, 
including  booms,  conical  surfaces,  etc.,  has  been  performed  only  to  the 
e.xtent  of  checking  for  reasonable  beh.avior  of  results  and  observing  num- 
erous '.’OMipi.'t  I'!’  ge''e’'ated  ipovic's  ('f  p:'i'*  i'-i'.'  t ra  i ect  I'r  i es  . |)irect  com- 
parisons I'f  results  uitti  another  code  shoiiKl  [Mobatily  be  made  here. 

Some  investigations  conducted  using  ABORC,  in  aduition  to  those 
given  in  .Appendices  B ,)iui  C,  an*  di^crilH-d  in  Re  ferenc  i-s  ')  through  11  . 
These  present  examples  of  various  cedi-  cajia!)  i 1 i t i es  such  as  the  modeling 
of  t iii’e  dejieiident  spectra,  re  entrant  geoiin't  r i i‘s  , and  cu  rren  t - i n j ec  t i on 
simulation  tests  of  sat  e 1 1 i ti'S  . 

, r.  Winaas  et  al,,  "Sensitivity’  ol'  SCf.MI’  Ri'sponse  to  Input  I'aram 
eters,"  !M.l,  Trans.  Nuc)  .'■'>c  i . .N'S-.’.',  .Vo.  (i.  lU-cr-mlHT  IP'S,  p.  ..'■■.(id. 

-’a,  ,1.  Wo, ids  aiul  R.  II.  Stahl,  "SCI  .'‘li’  Res]H)ns(,'  Sensiti\'it\  to  Ixplod- 
ing  Wire  Smirce  i’lioton  Output  \'a  r i ,i  t i oiis  . " I.Niid.  R’l  Sill  IU)2,  Sep  t emhe  i' 

Ls,  ip;:.. 

**’a,  .1.  Wouds  .iiui  !.  .\ . Pelmer,  "SK\’net  l'rov;i,im:  tin  lent  Inject  ion 
I’red  i c I ions  ■' So  1 lit  1 OH'i  to  .’'la  xwe  1 1 ' s 1 ijuaf  ions  , " li.VA  .■i!,,sp  1 . Iebrii,ir\  7, 

**l.  .\ . Pelmer,  "Response  ot'  a 1’ rech.i  rge  S.iti'llite  to  .i  Potent  i.il 
Photon  .’'O'.irce,"  I \TI  1.  Rl  Sill  P17.  N'ox’iaiiBer  71,  1 . 


Table  2 

SUMMARY  OF  ABORf  CODE  CAUART LITIES,  PHYSICS,  AND  MODE!, INC 
Geometrical  Configu rat i on 

Arbitrary  rotationally  symmetric  bodies  inside  arbitrary  rotationally 
symmetric  enclosures 

Arbitrary  energy,  angular,  and  axisymmetric  spatial  distributions  of 
emission  current  with  arbitiary  time  histories  can  be  specified;  ail 
disti ibur ions  can  be  time-dependent 

Two-dimensional,  axisymmetric  electron  orbits 

Two-dimensional,  axisymmetric  electric  and  magnetic  fields 

Limited  External  circuits 

Physical  Piieuomena  In'-luded 

Photocurrent  emissions  are  specified  to  the  code;  can  be  time- 
dependent  in  all  distributions 

Electric  and  magnetic  fields  and  current  densities  are  calculated 
Both  perfect  and  imperfect  conductors  can  be  specified 
Limited  dielectric  caj)abirity 

Fields  act  on  particles  of  charge  via  relativistic  equations  of  motion 

Time  dependence  treated  by  updating  Maxwell's  equations  each  time 
step  for  new  particle  positions 

Vacuum  only 

Method  of  So  1 ut  ion 

limission  current  broKen  int(j  particles  of  charge  with,  velocities; 
random  emission  in  all  variables 

Particles  of  ciiargc  are  converted  into  current  densities,  on  the  spa- 
tial grid 

tiurreiit  densities  are  list'd  to  obtain  E-  and  il-fields  b>'  I'initc- 
differencing  and  updating  Maxwell's  equations;  variable  ;:..;:'ne  sizes 

lilectric  ant!  magnetic  fields  act  hack  on  The  pai'ticles  via  the 
equations  of  motion,  altering  tlveir  tra  Jeetoj'i  es 
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DESCRIPTION  OF  THE  PHYSICS  AND  MODELING 
EMPLOYED  IN  THE  ABORC  CODE 


A-1.  INTRODUCTION 

Tliis  appendix  contains  a mathematical  description  of  the  ABORC  code, 
which  has  been  developed  to  nun;erjcally  evaluate  electron  motion  and  elec- 
tromagnetic field  generation  in  the  vicinity  of  arbitrary  bodies  of  revo- 
lution. The  problem  characteristics  are  taken  to  be  rotational ly  symmetric, 
reducing  the  problem  to  two  dimensions.  The  electromagnetic  fields  are  cal- 
culated from  Maxwell's  equations  and  are  used  to  influence  the  electron 
mot  ion . 

Electron  emission  from  the  surfaces  is  the  source  termCs)  for  driving 
the  problem.  Currents  and  fields  are  the  quantities  which  result  from  the 
calculation.  TliC  electron  emission  must  be  specified  in  space  and  time. 

For  example,  in  the  case  where  the  emission  is  due  to  photon  interaction 
with  materials,  the  photon  energy  and  time  spcctruni  determines  the  emis- 
sion charactei ist ics  of  the  electrons.  In  cases  sucli  as  this,  a separate 
electron  emission  code  must  be  used  to  get  the  electrcn  siiectnim  from  the 
photon  spectium.  (Particles  arc  used  to  represent  large  numbers  of  elec- 
trons.) The  quantities  calculated  directly  include  the  currents  in  the 
regions  between  the  conducting  surfaces,  as  well  as  the  electric  and  mag 
netic  fields  in  those  regions.  Ajijilying  Maxwell's  equations  properly  at 
the  boundaries  gives  surface  curre-its  aiul  cliarge  densities. 

The  remainder  of  this  appendix  outlines  the  methods  of  emitting  the 
jiariicles  of  charge,  calculating  [>article  motion,  converting  that  motion 
to  currents,  and  solving  the  field  equations.  The  field  equations  arc 
written  mi  a gene  ral  i ceil  coordinate  system,  so  the  conversion  from  the 
polar  coordinates  to  tlie  general  system  is  also  discussed.  i’olar  coordi- 
nates are  used  for  The  particle  motion. 
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I \-2.  ABORC  PARTICLE  EMITTER 

I 

I The  particle  emission  sections  of  ABORC  have  been  designed  to  allovv 

I for  relatively  easy  and  flexible  descriptions  of  complex  emission  charac- 

f teristics  and  to  simulate  "real-world"  electron  emission  with  a minimum 

I number  of  particles.  This  has  been  achieved  via  the  use  of  a free-form 

1' 

I input  processor  which  operates  upon  arbitrarily  specified  emission  descrip- 

I; 

I tions,  and  by  implementation  of  general  Monte  Carlo  methods  to  the  entire 

S range  of  emission  characteristics  — i.e.,  intensity,  energy,  angle,  space, 

and  time.  We  describe  herein  the  features  of  the  emission  specification 

; method  and  provide  examples  of  the  use  of  the  emitter.  A description  is 

included  of  the  operation  of  the  emitter  itself,  the  methods  employed, 
and  the  benefits  derived. 

A- 2 . 1 ABORC  Emission  Input  Description 

The  problem  of  describing  to  the  computer  code  exactly  what  emission 
is  to  be  simulated  is  twofold:  (1)  description  of  the  particle  character- 
istics and  intensity,  and  (2)  the  position  in  the  computer  model  where 
this  emission  is  to  take  place.  The  following  describes  the  essential 
features  which  interface  the  analyst's  desires  with  the  ABORC  internal 
arrays,  flags,  equations,  etc. 

F:  SION  CHARACTERISTICS 

The  above  card  simply  informs  the  code  that  emission  information 
fol lows . 

EMISSION  INTENSMY  n --  POINT  PAIR  TABLi 

The  emissior!  intensity  ^'ormat  is  used  to  specify  the  time  history  for 
the  emission  pulse.s;  n is  the  intensity  number,  of  which  there  can  t)e  uji 

7 

to  20,  anti  POINT  PAIR  TABLE  is  the  time  history  in  anip/n;“  versus  time  in 
sec.  El  may  he  used  to  abbre>'iatc  "emission  intensit)'."  i'tu'  exnirijsle,  a 
2()-n:'.ec  lAVHM  trCmgular  pulse  wuth  a peak  of  0.1  amp/iii*  wt.iuld  be  input: 

b = 0,0,20[>9.  . 1 ,2  01-9,0 
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I Ml  00  ION  INIINSIIY 
(or  I r:o 


The  table  may  contain  up  to  50  point  pairs,  allowing  high  resolution 
to  fairly  complex  emission  time  histories,  and  may  be  continued  on  as  many 
cards  as  desired. 

ENERGV  DISTRIBUTION  n (BINS  = i)  - POINT  PAIR  TABLE 

Similarlv,  the  energy  distribution  card  allows  for  the  arbitrary  spe- 
cification of  the  emission  energy  spectrum;  n is  the  distribution  number, 
and  again,  up  to  20  are  allowed  at  present.  FD  may  be  used  to  abbreviate 
"energy  distribution."  (BINS  = 1)  is  optional.  i specifies  that  there 
shall  be  i particles  emitted  per  particle  time  step  per  emission  zone  for 
the  distribution.  If  the  (BINS  = i)  input  is  omitted,  the  default  value 
of  i is  4.  The  POINT  PAIR  TABLE  is  as  before,  with  relative  frequency 
fnumber/uni t energy  versus  energy  (eV)]  as  the  parameters.  For  example, 
the  spectrum  of  Figure  A-1  might  be  input: 

ENER(3Y  DISTRIBUTION  7 * 4K,  0,  5K, . 4,6K, . 8 , 7K,  1 , 

9K,  .8,  12K,  .4, 15K,  . I5,20K,0 

Mono-energetic,  emission  is  specified  simply  by  EDn  = energy. 

Note  the  use  of  the  "K,"  which  is  interpreted  as  10'^.  Various  other 
letters  may  be  used  with  floating  point  numbers  in  the  emission  section; 
they  are  listed  in  Table  A-1.  The  example  spectrum  is  given  normalized 
to  a peak  of  1;  however,  the  emitter  normalizes  all  distribution  functions, 
so  the  emission  current  density  has  the  value  specified  by  the  "TIMES"  fac- 
tor described  bolo^^  under  "EMISSION  ZONES,"  and  multiplied  by  the  "FLUENC" 
factor  discussed  in  Ajipendix  E.  Desired  relative  height  for  all  distri- 
butions is  all  the  user  need  consider. 

ANGLE  DISTRIBUTION  n - POINT  PAIR  TABLl 

The  angle  distribution  card  is  shown  aboi'o.  AD  may  be  used  to  abbre- 
viate "angle  distribution."  The  parameters  for  the  point  pair  table  in 
this  instance  arc  the  ;iumbcr  of  jrarticles  per  unit  angle  versus  angle  in 
radians.  For  example,  the  distribution  of  Figure  A-2  is  given  by; 

ANGLE  niS'TRIBUTION  5 -■  0 , 0 , . 6 , . 7 b , 1 . 1 , . 45  , 1 . b / , 0 

(or  AD3) 
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Figure  A-1.  Sample  energy  spectrum  input  to  ABORC 


Table  A-1 

FLOATIN  ; POINT  MULTIPLIERS  FOR  ABBREVIATED  ABORC  INPUTS 


Letter  Factor 

Power  of  10 

P 

-12 

N 

-9 

U 

-6 

K 

+ 3 

M 

+ 6 

The  letter  factor  can  be  used  to  replace  the  power  of 
10  desired  on  input  cards. 
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Figure  A-2.  Sample  angle  distribution  input  to  ABORC 

Polar  angles  are  employed  in  the  coding  for  emission  specification.  The 
angle  on  the  given  sample  is  the  angle  from  the  emitting  surface  normal. 
The  azimuthal  angle  distribution  is  specified  uniform  from  0 to  2tt  for 
two-dimensional  calculations.  Straight-out  emission  is  specified  simply 
by  ADn  = 0. 

EMISSION  ZONES 

The  emission  zones  card  simjjly  specifies  that  emission  zone  informa- 
tion follows.  The  emission  zones,  then,  are  indicated  by: 

(2^,ry  TO  ^Z2'’"2'  ^ k,ADQ3  I 

Zj,  I'j  , z^,  give,  of  course,  tlie  coordinates  of  the  zone.  i is 
tl\e  emission  intensity  time  history  number  reference  to  be  used  for  the 
zone,  delayed  l^y  x seconds  and  multiplied  by  y.  The  DELAY  and  1 IMPS  tac- 
tors  are  optional.  j is  the  energy  distribution  number  to  be  used  for 
ttie  zone,  k is  the  angle  distribution  to  be  used  for  dcterniining  the 
direction  of  cacti  jiarticle  relative  to  the  surface  normal,  and  I is  the 
angle  d i st  ri  tiut  ion  refei'ence  for  the  emission  electron  direction  azimutlial 


angle  about  the  surface  normal.  Particles  are  always  emitted  to  the  left 
of  the  line  segment  running  from  (Zj,rp  to  Therefore,  care  must 

be  taken  to  specify  zones  in  proper  order. 

As  an  example,  suppose  we  desire  an  emission  zone  position  on  the  top 
of  a cylinder,  from  the  axis  to  0.2-m  radius  and  at  vertical  position  0.5 
m.  We  shall  reference  the  previous  distributions  and  input: 

(.5,0)  TO  (.5,. 2)  INT  5 TIMES  10,  ED7,  AD3,  ADQ3  2 

vjhere  we  have  multipled  the  current  density  by  10.  We  may  also  want  to 
emit  from  the  side  of  the  cylinder  (e.g.,  radius  = 0.2  m)  but  with  reduced 
intensity  and  delayed  by  a couple  of  nanoseconds: 

TO  (.3,. 2)  INT  5 DELAY  2E-9  TIMES  .1,  ED7,  AD3,  ADQ3  2 

where  z^  and  r^^  for  this  emission  zone  are  obtained  from  the  previous  zone's 
z^  and  r^. 

END  OF  INPUT 

terminates  the  emission  input  processing, 

A- 2 . 2 ABORC  Emitter  Operation 

The  emitter  operation  is  outlined  here.  Generally,  the  idea  has  been 
to  randomize  the  entire  emission  process  instead  of  emitting  from  discrete 
points  in  space,  angle,  energy,  and  time,  as  has  been  done  in  the  past. 

This  randomization  reduces  the  systematic  excitation  of  high-frequency 
modes,  and  represents  physically  continuous  distributions  more  accurately 
than  discrete  points  for  the  same  number  of  particles.  This  treatment  of 
emission  currents  is  particularly  beneficial  under  high  SCL  conditions, 
where  a small  fraction  of  the  emitted  charge  may  be  causing  dominant 
response . 

The  emitter,  then,  based  upon  the  amount  of  charge  to  be  emitted  at 
a given  time  step,  sets  up  particles  with  randomly  determined  energies 
from  the  given  energy  spectrum,  randomly  determined  angles  (velocity  com- 
ponentsj  consistent  with  the  given  angular  distribution,  randomly  deter- 
mined position  in  the  emission  zone  (taken  to  be  uniform  over  the  zone 
area),  and  randomly  determined  emission  time  (taken  to  be  uniform  ove>' 
tlie  particle  time  step). 
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The  energy  characteristics  of  the  emitted  charges  are  taken  to  be  cru- 
cial to  effective  simulation.  For  coding  purposes,  this  translates  to 
allowing  very  precise  definition  of  the  desired  spectrum  (discussed  in  the 
input  section)  and  effectively  sampling  values  from  this  spectrum  without 
generation  of  too  many  particles  (which  could  result  in  excessive  running 
times).  We  achieve  these  objectives  via  the  use  of  the  "stratified  sam- 
pling" scheme.*  Given  a spectrum  such  as  that  presented  in  Figure  A-1, 
we  transform  (via  integration/normalization)  to  a distribution  function 
of  the  form  shown  in  Figure  A-3. 
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A comparison  of  "stratified"  sampling  with  "straight”  samjiiuig  - i.e. 
sampling  from  the  entire  distribution  insteac  of  I'rom  ranges  within  tiic 
distribution  — will  iliu:;trate  the  benefits  of  stratification.  Straiglit 
Monte  Carlo  sampling  applied  to  this  distribution  (i.e.,  gcin'rat ion  ot  a 

. liammersby  and  I',  llandseomb , Monte  Carlo  Methods,  .! . Wiley  (|  dons, 
New  York  ( 196B) . 
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random  variate  between  0 and  1 and  interpolating  the  function  for  the 
corresponding  energy)  yields  a standard  sampling  error  s: 

n 


i=l 


where 


p = mean, 

e^  = i^^  sampled  energy, 
n = number  of  samples  taken. 

The  sampling  error,  then,  goes  inversely  with  the  square  root  of  the  number 
of  samples,  or  four  times  as  many  samples  must  be  taken  to  increase  the 
accuracy  of  s by  a factor  of  2. 

The  use  of  stratified  sampling  (i.e.>  dividing  the  0--to-l  interval  into 
k strata  and  sampling  from  within  each  sti'atum)  gives  a standard  error: 


(Aa.)' 


z 


j = l ■'  ■’  i = l 


(e.  . 
ij 


h) 


2 


where 


k - number  of  strata, 

= random  variate  range  for  the  interval, 
iij  = number  of  samples  taken  in  the  interval. 


c.  . 
].) 


energies  selected. 

't  is  seen  that  the  error  goes  inversely  with  the  number  of  s;uiiples  taken; 


i.e.,  doubling  tiie  number  of  samples  roughly  halves  the  standard  error. 
This  translates  to  better  .statistics  for  fewer  particles. 

The  directional  cnaracteri st ics  of  the  emitted  charge  are  important 
also.  lixperieiu'e  indicates  that  rhe  angle  di  str  itnit  ion  is  not  as  criti  - 
cally iiajjortant  as  tiie  pr.rper  ei’orgy  a i st  r i liut  ion  , particular’}'  ;\t  higher 
fluence  conditions.  Wc  allow,  then,  the  suec i t i cat  i on  of  anguiar  charac - 
teristics  and  apply  standard  Monte  (,'arlo  saiii|iling  to  these  d i .st  r i but  i on  s 
instead  of  the  more  sophisticated  sti’atified  variety. 

With  rogar'd  to  the  sj'ace  and  time  characteristics,  we  take  tine  emis- 
sion to  he  uniform  o\'.e;'  the  emission  eone  area  aiul  ova'i'  the  emission  time 


;ss 


Step.  Since  the  intent  is  to  simply  spray  the  emitted  charge  from  the  zone 
at  random  positions  in  the  space  interval  and  at  random  times  in  the  time 
step  interval,  we  find  that  the  standard  Monte  Carlo  is  adequate  here  also. 

In  summary,  the  emitter  calculates,  from  the  specified  current  den- 
sity time  history,  the  charge  to  be  emitted  at  each  emission  time  step  in 
each  emission  zone.  The  energy  characteristics  are  randomly  and  accurately 
selected  using  the  stratified  energy  distribution  sampling.  The  particles 
are  then  assigned  directional  characteristics,  emission  time,  and  position 
via  standard  Monte  Carlo  sampling. 

An  outstanding  capability  of  this  particle  emitter,  aside  from  the  ease 
of  inputs  and  specification  of  complicated  characteristics,  is  its  ability 
to  model  time- dependent  energy  spectra.  This  capability,  necessary  for 
complete  modeling  of  photo-electron  emission  results  obtained  from  photon- 
generating machines,  is  available  by  using  the  existing  coding  and  simply 
specifying  different  energy  distributions  with  different  time  histories 
for  emission  from  a common  spatial  zone.  As  many  as  20  energy  distribu- 
tions and  20  time  histories  can  be  specified  to  model  the  complete  time- 
dependent  spectrum,  if  necessary. 

A- 3.  ABORC  PARTICLE- PUSHER 

The  finite  "particles"  of  charge  emerging  from  the  emitter  portion 
of  the  code,  discussed  in  the  previous  section,  are  picked  up  and  prop- 
agated through  the  spatial  zones  with  forces  acting  on  tliem  due  to  the 
E-  and  li-fields  in  the  space.  The  jnirticlc  trajectories  are  determined 
by  the  equations  of  motion,  whicti  are  updated  each  pai'ticle  time  step. 

The  old  and  new  particle  iiositions  are  then  employed  to  generate  currents. 
The?  currents  are  used  to  update  the  E-  and  !l-fields.  A descrijition  of  the 
determination  of  particle  trajectories  in  ABORC  is  given  here. 

The  particle  has  rest  mass  in  and  chart'!-  ((,  and  is  in  an  electric  field 
i;  - iE.,,Ej,,0)  and  magnetic  field  B = (0,1), Bp).  The  eipi.jtion  of  motion  can 
lie  numerically  integrated  in  two  steps; 

1.  Eind:  the  now  magnitude  of  the  -yelocily  v from  the  energy 
equal  i oi.s  . 

2,  I'ind  the  direction  of  t!ie  velocity  frum  momentum  equations. 


The  equations  were  originally  updated  by  differencing  the  momentum 
equations  in  cylindrical  coordinates.  This  had  the  effect  of  introducing 
a singularity  at  the  axis  of  the  coordinate  system.  To  eliminate  this 
singularity,  the  particle  motion  is  now  described  in  Cartesian  coordinates. 
This  was  a simple  and  expeditious  means  of  solving  the  problem. 

The  force  is  calculated  from  the  nearest  field  points  of  the  grid  in 
the  following  manner.  The  magnetic  field  is  equal  to  the  field  at  the 
nearest  point.  Since  in  SGEMP  problems  the  magnetic  force  is  small,  this 
is  sufficiently  accurate.  For  the  electric  fields  — say,  in  the  Z direc- 
tion — the  two  nearest  fields  in  that  direction  are  found  to  have  the  same 
coordinates  in  that  direction.  These  two  are  averaged  to  get  the  field 
acting  on  the  particle.  This  has  the  advantage  of  giving  forces  on  emitted 
particles  which  are  given  by  the  normal  forces  just  outside  the  emission 
surface  even  though  in  free  space  there  are  self-forces  on  the  particles. 
(This  latter  fact  may  be  significant  for  bodies  with  booms,  where  a long 
distance  must  be  traversed  by  the  electrons.) 

For  motion  of  the  particles,  some  confusion  exists  due  to  the  time 
centering  of  the  problem.  A primary  cause  of  this  is  the  fact  that,  as  a 
particle  is  injected  into  the  problem,  the  first  time  step  is  important 
since  violent  motion  occurs  in  this  step  for  low-particle  energies.  We 
use  a random  time  step  for  the  injection  so  that  low-energy  ixirticles  do 
not  systematically  leave  the  |)roblem  in  one  step;  tlierefotu',  the  centerijig 
is  not  strictly  defined.  Still,  we  liave  found  the  following  sctieme 


adequate . 

Let  the  superscri|)t  " + " denote  new  v.ilues  and  \ the  posiricn  and 
tlie  velocity  of  ttie  particle.  i'he  |)aiticle  is  treated  as  hi.‘ing  in  tlir(.‘e 
dimensions  aiul  rotateil  about  the  .ixis  of  s>’mmetr\'  back  to  .j  coordinate 
system,  where  one  of  t lie  coordinates  i •;  ceru  to  save  storage. 

Using  the  not  at  i cm 


i 


where  c is  the  cclocit>'  of  light, 


v'*'  = V + 3^  (^  + V X ^) 

my 


+ qAt  ->  ,->■+ 

Y = Y + E*iv  + vl 


E* (v  + v)  . 


Making  one  iteration, 


= y-x  + 331  [H  ^ o.5(r  -H 1)  X Sj  . 

Y mY 

y'^'^  = Y + -3^  [E-(v'^‘^  + v)]  , 


-'■+  ->  ->++  , 

X ■-  X + V At, 


where  q and  m are  the  charge  and  mass  of  the  electron  and  At  is  the  time 
step.  This  oneration  is  followed  by  the  abovementi oned  rotation  which 
involves  taking  a square  root  and  is,  therefore,  not  ideal. 


A-4.  DESCRIPTION  OF  CllAROED-PARTICLE  MOTION  BY  CURRENTS 

In  the  numerical  solution  of  Maxwell's  equations,  currents  are  eval- 
u.'itetl  at  discrete  points  in  sjiace  and  time.  I’articles,  representing  elec- 
trons, move  in  a generalized  coordinate'  system  in  time.  An  interpolation 
scln.'ine  must  lie  constructed  to  translate  the  continuous  motion  of  the  ])ar- 
ticles  into  a set  of  discrete  currents,  The  scheme  chosen  is  such  that 
the  residual  charge,  as  calculated  from  the  time  integral  oi'  the  diver- 
gence of  the  current  is  'ero  after  a particle  lias  passeti  into  and  tlien 
out  of  .1  region.  (This  does  not  lioUi  fo'  boundary  zones  since  currents 
oi.'tside  tile  region  of  interest  are  not  eonsidered.) 

In  the  present  zoning,  currents  aie  evaluated  at  .'oin'  centers  in 
their  own  direction  and  zone  bouiid.i  r i es  in  perpendicular  directions. 

!'l,  ? currcnit  is  . inie- centereU . .Accordingly,  as  forci'd  lo'  tlie  continuity 
equation,  t lie  chargi'  dens'.ties  .ire  at  zone  boundaries  in  space  and  time. 

Ihe  present  obiective  is,  then,  to  take  .i  particle  from  point 


.!  a a .1 
(Mj.q_,  .q-,t  ) 


to  the  point 

^ b b b .b. 

Id  9 

where  t = t + At.  That  is,  the  currents  representing  such  a translation 
must  be  generate^. 

Since  the  grid  spacing  is  uniform  in  the  q's,  linear  interpolation 
can  be  performed  in  the  transverse  coordinates.  This  can  be  seen  from  the 
formulation 
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r b 


- 'in 


At 


Q.,  (i+% 


. j ,k)j 


Q|  Aq^  (i^kJ/K)  Aq,^  (i  + ’2,j,k)  Aq^ 


for  point  (q.,q,  q^J  = [ (i+'2)Aqj , jAq^  ,kAq^]  , particle  of  charge  . Hc’"e, 

J is  the  current  in  the  q,  direction. 

I ». 

The  first  bracket  contairis  the  real-space  particle  velocity  in  the  q, 
direction.  ftie  second  bracket  contains  the  real -space  volume  element.  The 
integral,  represents  the  t ime-as'erage  fractional  distance  from  the  point  of 
interest  in  the  transverse  directions.  It  can  be  seen  that  th.e  particle  is 
treated  as  a volume  element  in  q space  of  dimcnsi..ns  Aq^,Aq^,Aq,,  and  the 
integral  represents  the  tra.'isverse  area  overlapping  the  cone  of  interest. 

The  .lime  limits  arc  such  that  they  are  within  the  range  of  interest, 
t^\t‘\  and  the  particle  is  within  one  zone  of  the  point  of  interest  in  its 
linear  traversal: 


a 


, , b 

(Ml 


q‘i 
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from  point  a to  point  b. 

In  practice,  a particle  is  tracked  from  point  a to  point  b in  inter- 
vals of  zone  crossings  so  that  the  appropriate  values  of  i,  j,  and  k are 
easily  established.  The  treatment  for  <and  is  related  Co  that  of 
by  the  same  synunetry  as  mentioned  in  the  discussion  of  the  field  calcula- 
tion below. 


A-5.  SOLUTION  OF  THE  FIELD  EQUATIONS 


The  numerical  solution  of  Maxwell's  equations  is  straightforward. 
The  initial  conditions  are  that  all  fields  and  charge  densities  are  zero 
The  currents,  as  calculated  from  the  motion  of  the  charged  particles 
injected  into  the  region  (.^ee  below),  are  the  quantities  which  drive  the 
time  evolution  of  the  fields.  Thus,  for  a medium  with  the  permittivity 
(Cq)  and  permeability  (y^)  of  free  space.  Maxwell's  equations  reduce  to 

=0 1!  = ^ ” 


and 


with  tlie  initial  conditions 

^ = 0 at  t = 0 , 

M = 0 at  t - 0 , 

^ X n - 0 where  n is  the  normal  to  the  bounding  surface. 

Putting  these  equations  imo  numerical  fo"m  for  solution  is  also 
straightforward,  and  may  be  done  directly  in  two  dimensions.  However, 
by  jrerforming  the  task  in  three  dimensions,  one  is  forced  into  a symmetry 
which  is  very  convenient  and  not  obvious  in  the  two-dimensional  case. 

To  complete  the  geometrical  generality  (and  permit  a simple  method 
of  varying  zone  spacing),  tlie  sjjace  under  consideration  is  taken  to  he 
metrized  by  the  generalized  orthogonal  coordinates  (q^.q^.q.^)  where  the 
order  is  ‘uch  tnat  the  coordinate  system  is  right  handed.  Using  the 


40 


•k 

notation  of  Margenau  and  Murphy,  displacements  in  real  space,  ds,  may  be 
related  to  displacements  in  q space  by  the  functions  Q,  where 

ds  = Q.  dq^  . 

Using  this  definition  and  the  definition  of  the  curl,  Vx,  as  the  path 
integral  of  a quantity  in  the  left-hand  direct;'Ofi  (Cauchy  rule)  around 
a closed  path  divided  by  the  area  of  the  enclosed  surface,  Maxwell's 
equations  in  finite-difference  form  can  be  represented  by 


E,  (1+^,1, h)  = E^  ('+'2,j,k) — J,  (i+-S,j,k) 

'0 


At 


I 

2:^  Q,  (i^^iTT.kj  Aq^  Q,  (i^j.k^  Aq^  | 


x'.iH,  ' M + ‘2,j+-2.k)  (i+’2,i+%,k)  Aq^  - ( 1 + -^, j --d, k)  Q.^  (i+b,j-%,k)  Aq 


(i+‘'2,  j ,k+b)  (i+''2,j  ,k+J2)  Aq,,  - ^(i+-i,j  ,k-h)  Q.,  (i+''2,  j .k-^'i)  Aq 


and 


(i,  j+'^/K+Js,  = {i,]+hMh) 


At 


n+J 


Q,  (k,  j+'a,k+'5)  Aq^  Q,  ( i , j +'-2, k + ’i | Aq^ 


(i  , j+'2,k  + r)  Q.j.  (i  . j + 1 ^k  + 'i)  Aq^^  - li'^'' ^ {i,j,k  + 'f)  Q,,  (i.j.k^'i)  Aq^ 

e”'’^  (i,j+2,k  + l)  (.  i , j+'2,k  + l ) Aq^  - {i,j+'2,k)  Q,,  (i.j+'i.k)  Aq,, 


vvi,  th  the  continuity  equation  (not  essfuitial  to  determining  tiie  iield^,  ljut 
of  interest  for  itself  and  for  converting  particle  morion  into  equivalent 
currents)  rejn-esented  by 


il.  Margenau  and  C.  M.  Murphy,  Tlu'  M.it  hema T i c s of  l’Ii>'sic:^  and  Clu  iii- 
istrv’,  Princeton,  1'.  van  Nostrand  Company,  Inc.  (1956). 
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n+l  , . ...  11  , . , . 

p (1,3  JO  =-■  P (1.  J ,k) 


At 


n+% 


(; 


C^  (i,j,k)  Aq^  q^  ri.j>k)  Aq^  Q^  fi,j,k)  Aq^ 

j|j“'''''  (i+5'2,j,k)  (i+^i.j.kl  Aq^  Q^  (i+V.,j,k)  Aq^ 

- j”’’^  (i-H.j.k)  (i-li,j,k3  Aq^  Q^  (i-'j.j.k)  Aq^ 

(i,j+*1,k)  Q^  (i,3+  %k)  Aq^  Q^  [i+^2,j,k)  Aq^^ 

- (i,j-i5,k)  Q^  (i,j“?2,k)  Aq^  Q,  (i.j-’j.k)  Aq^j 


(i,j,k+52j  Qj  (i,j.k+%)  Aq,  Q^  (i , j , k-i-’^)  Aq^ 

- (i,j,k-}i)  Q^  (i,j.k-!-i)  Aq^  (i,j,k-’i.)  Aq^  | . 

J ) 

The  notation  used  in  these  equations  requires  some  comment.  The  super 
script  refers  to  the  time  step.  Thus,  some  quantities  are  centered  in  time 
and  some  are  at  boundaries  in  time.  Corresponding  to  this,  there  are  two 

time  steps,  the  time  step  connecting  quantities  centered  in  time  At'\,  ind 

n ' 

the  time  step  connecting  quantities  at  boundaries  in  time,  At 


Subscripts  refer  to  directions  in  the  generalized  coordinate  space. 
Thus,  is  the  component  of  the  ciectric  field  along  the  direction  of  a 
displacement  in  space  given  by  a displacement  in  q,  at  tiie  spatial  point 

i 

in  que.'-tion.  Quantities  in  parentheses  refer  to  the  iiosition  in  space. 
Thus,  is  evaluated  at  the  point  in  space  determined  by  the 

coord  in  at  es 


M]  = « + ‘liCMlN)  ’ 

q,  = P Aq.^  + ‘JlMIN)  ’ 

where  tiie  liiiniinum  value  of  the  coordinate  is  specified  for  convenience, 
allowing  the  -.'.patial  boundaries  to  be  other  tiian  zero  in  the  q space. 


It  will  be  noted  thnt  the  grid  spacing  in  q space  is  uniform.  Fur- 
ther, if  one  of  the  Q's  is  zero  at  a point  of  interest,  the  procedure 
fails.  In  fact,  at  such  points  the  coordinate  system  does  not  metrize 
real  .space.  The  failure  is  that  many  points  in  q space  correspond  to 
one  point  in  real  space.  Such  cases  must  be  treated  specially. 

The  equations  for  the  other  components  or  the  electric  and  magnetic 
fields  are  obtained  by  cyclically  permuting  the  integer  subscripts  and 
the  corresponding  coordinates  a,  6,  and  y- 
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ABSTRACT 


General  scaling  laws  have  been  previously  derived  and  used  for 
numerous  applications  in  electromagnetic  theory.  The  purpose  of  this 
appendix  is  to  report  investigations  of  the  scaling  laws  specifically 
for  the  nonlinear  system-generated  electromagnetic  pulse  (SGEMP)  problem, 
and  to  use  the  scaling  laws  to  present  parametric  SGEMP  calculations 
applicable  to  a wide  range  of  pertinent  excitation  parameters  including 
pulse  width,  fluence,  energy  spectrum,  and  object  dimensions. 

According  to  the  scaling  laws,  if  the  pulse  time  history  and  object 
dimensions  are  scaled  by  a factor  a,  the  incident  photon  fluence  is 
scaled  as  1/a,  and  the  emitted  electron  energy  distribution  is  unchanged, 
the  resulting  electromagnetic  response  of  a perfectly  conducting  body 
will  scale  as  shown  in  Table  B-1. 

Results  of  specific  SGEMP  calculations  are  presented,  using  the 
ABORC  code  in  such  a format  as  to  make  possible  structural  response 
estimates  for  a wide  variety  of  conditions,  Interesting  trends  in  the 
response  of  simple  objects  as  a function  of  excitation  parameters  are 
observed  and  discussed. 

B-1.  INTRODUCTION 

General  scaling  laws  have  been  previously  derived  (Refs.  1,2)  and 
used  for  numerous  applications  in  electromagnetic  theory.  This  appendix 
investigates  the  scaling  laws  specifically  for  the  nonlinear  system- 
generated electromagnetic  pulse  (.EGEMP)  problem  and,  using  the  scaling 
laws,  presents  parametric  calculations  applicable  to  a wide  range  of 

iGeorge  Sinclair,  "Theory  of  Models  of  Electromagnetic  Systems, 

Proc.  IRE,  Vol.  36  (1948),  p.  1364-1370. 

2t.  N.  Delmer  et  al . , "SGEMP  Phenomenology  and  Computer  Code 
Development,"  DNA  3653F,  November  11,  1974. 

3e . P.  Wenaas,  S.  Rogers,  and  A.  J.  Woods,  "Sensitivity  of  SGEMP 
Response  to  Input  Parameters,"  IEEE  Trans.  Nucl.  Sci . NS-22,  December  1975. 
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pertinent  excitation  parameters  (Ref.  5)  including  pulse  width,  fluence, 
energy  spectrum,  and  object  dimensions.  Interesting  trends  in  the  response 
of  a simple  object  to  variations  in  the  excitation  parameters  are  noted 
and  discussed. 


B-2.  SCALING  LAWS 

The  quantities  most  readily  scaled  are  time  t,  object  dimension  r, 
incident  photon  fluence  $,  and  electron  energy  or  velocity  v.  To  determine 
how  electromagnetic  quantities  of  interest  scale  with  t,  r,  $,  and  v,  one 
need  only  substitute  the  scaled  quantities  shown  below  in  the  applicable 
equations  describing  the  electromagnetic  response. 
t'  = t/r  = r<t> 

r"  = r/R  v'  = Bv 

For  the  linear  regime  where  the  electron  motion  is  unaffected  by  the 
electric  and  magnetic  field  forces,  the  set  of  Maxwell's  equations  is 
sufficit'  c to  describe  the  electromagnetic  response  of  a perfectly 
conducting  body.  Substituting  the  scaled  quantities  into  Maxwell's 
equations  results  in  the  requirement  that  the  nondimensional  quantities 
T and  R be  equal  and  that  the  quantity  6 be  unity.  The  resulting  scaled 
electromagnetic  quantities  for  the  linear  regime  are  shown  in  Table  B-1, 
where  R = r A a . 

If  the  electron  trajectories  are  affected  by  the  electric  or  magnetic 
fields  (nonlinear  regime),  then  Newton's  first  law  describing  electron 
motion  under  the  influence  of  electric  and  magnetic  field  forces  must  also 
be  considered.  The  additional  requirement  resulting  from  the  scaling 
substitution  into  Newton's  law  is  that  the  nondimensional  parameters 
r and  a must  be  equal.  The  resulting  scaled  electromagnetic  quantities 
applicable  to  the  nonlinear  regime  are  also  shown  in  Table  P-1. 
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Table  B-1 
SCALING  LAWS 


Quantities 

Scaled  Quantities 
(Linear  Regime 

Scaled  Quantities 
(Nonlinear  Regime) 

Current  density,  J 

J'  = aTJ 

.J'  - 

Charge  density. 

p'  = aVp 

. 2 

P = a p 

Electric  field,  E 

E'  = PE 

E'  r:  aE 

Magnetic  field,  H 

= FH 

H'  - aH 

Potential  difference. 

AV  AV'  = r/aAV 

AV'  = AV 

Current , I 

1'  = r/al 

r = I 

B-3.  RESULTS 

SGEMP  response  calculations  using  the  ABORC  computer  code  (Ref,  4) 
were  nerformed  parametrically  as  a function  of  the  various  excitation 
parameters.  Responses  resulting  from  a wide  range  of  excitation 
parameter  values  can  be  displayed  by  utilizing  the  scaling  laws. 

As  an  example,  calculations  ha^re  been  performed  for  the  response  of 
a right  circular  cylinder  of  radius  R exposed  to  a photon  pulse  of 
:^luence  having  a Planck  radiation  spectrum  characterized  by  an  energy  2. 
The  pulse  time  history  is  proportional  to  sin  [(-rit/Z) /t]  . Note  that 
becausi’  of  the  scaling  laws  one  need  not  select  specific  values  for 
all  the  various  excitation  parameters  to  display  the  results;  they  may 
be  displayed  parametrically  as  functions  of  t/2itR/c,  ifR,  and  2.  The 
nondimens ional  quantity  •r/(2rR/c)  is  simply  the  ratio  of  pulse  rise  time 
to  the  time  required  for  light  to  travel  around  the  object.  The  quantity 
$R  is  the  product  of  fluence  and  object  dimension  which  we  define  as  the 
"fluence-product,”  having  dimensions  of  cal/m  in  the  mks  system. 

For  purposes  of  this  summary,  the  surface  current  at  midpoint  on  tlie 
side  of  the  cylinder  is  the  quantity  chosen  to  characterize  SGEMP  response. 
Results  of  the  computer  calculations  are  shown  in  Figures  i<l  ami  I'.,’. 
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Te  B-1.  Surface  current  versus  fluence  product 
different  spectra  characterized  by  energy^in  keV 


In  Figure  B-1,  the  total  current  flowing  on  the  cylinder  across 
the  midpoint  is  shown  as  a function  of  the  flusnce  product  for  various 
photon  spectra.  The  following  observations  are  drawn  from  Figure  B-1. 

• For  characteristic  energy  in  the  range  of  1 to  100  keV,  the 

response  is  essentially  linear  for  f luence-products  (R'J’)  less 
_2 

than  10  cal/m.  In  this  range,  the  softest  spectrum  yields 
the  largest  emission  current  and,  therefore,  the  largest 
r 'sponse. 

• The  responses  became  nonlinear  with  increasing  fluence  products, 
and  the  harder  spectra  yield  higher  responses  because  the  more 
energetic  electrons  in  these  spectra  are  less  space-charge-limited. 

In  Figure  B-2,  the  total  current  is  shoi'n  as  a function  of  characteristic 
energy  for  various  fluence  levels.  The  following  observations  can  be  .lU’de. 

• .\t  the  lower  fluences  in  the  linear  regime,  the  response  falls  off 
as  1/S  which  corresponds  roughly  to  the  fall-o^-^  in  emission 
current  density  with  increasing  photon  energy. 

• At  higher  fluences,  there  is  a definite  peak  in  the  maximum 
response  as  a function  of  energy.  Thvi^,  for  a given  fluence, 

there  is  a defined  energy  range  whicl;  produces  a worst-case  resp>  ise 
The  previous  results  have  been  obtained  for  a constant  ratio  of  pulse  width 
to  object  dimension. 


B-4.  SUMMARY 


Scaling  laws  for  SGEMP  excitation  are  defined. 

SGEMP  calculations  have  been  performed  and  displayed  in  a compact 
format  using  the  scaling  laws.  Results  are  useful  for  making  order-of- 
magnitude  external  SGEMP  response  estimates  for  a wide  variety  of 
excitation  parameters . 
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SAIvIPLE  .ABORC  CALCUUTIONS 
SGEMP  GEOMETRY  EFFECTS 


ABSTRACT 


Self-consistent,  fully  dynamic  computer  calculations  were  performed 
using  the  recently  developed  arbitrary  body -of-revolution  code  ABORC  fr.r 
coiviiplex  geometries  in  SCSMP  environments  to  test  the  validity  of  simpli- 
fying geometry  assumptions  previously  made  in  the  solutions  of  these 
problems.  Assumptions  such  as  simple  geometry  representations  of  complex 
bodies  and  separability  of  inside  and  outside  problems  are  tested.  Effects 
of  gaps,  interior  electrical  paths,  and  booms  protruding  from  bodies  are 
discussed.  High  space-charge-limited  results  are  emphasized. 

The  response  of  simple  bodies  such  as  cylincers  and  spheres  was 
found  to  be  similar  to  move  complex  gecmetries  in  terms  of  peak  currents 
and  response  rise  times,  although  much  of  the  detail  in  terms  of  resonant 
frequencies,  internal  responses,  and  responses  around  reentrant  bodies 
is  lost.  The  inside  and  outside  responses  of  an  object  can  be  separated 
in  many  cases  even  though  both  solutions  by  themselves  are  highly  nonlinear 
and  the  leakage  currents  are  relatively  large.  The  external  response  of  a 
highly  segmented  body  is  similar  to  the  response  of  a smoothly  connected 
body  of  revolution  'Fhe  internal  response  can  be  considerably  different 
when  a conducting  path  exists  between  the  segments. 

C-1.  INTRODUCTION 


The  problems  associated  with  computing  the  SGIiMP  response  of  a complex 
structure  caa  oe  conveniently  divided  into  two  categories:  the  choice  of 
structural  configuration  that  best  reprasents  the  complex  object,  aiid  the 
choice  cf  various  input  excitation  parameters  foi  c given  geometry.  Problems 
assoc  ated  with  the  response  of  different  structural  configurations  are 
testei'  in  this  paper,  while  problems  associated  'nth  differences  in  excitation 
parameters  for  a given  configuration  are  treated  in  a companion  papei  (Ref.  1) . 


5.4 


The  structural  detail  that  can  be  treated  by  the  SGEMP  codes  has 
increased  from  simple  spheres  and  cylinders  (Refs.  2,3)  to  two-dimensional 
arbitrary  bodies  of  revolution  (Ref.  4) , and  techniques  are  being  developed 
to  treat  fully  sclf-ccnsistent  problems  in  three  dimensions.  The  increase 
in  modeling  detail  has  been  pursued  to  gain  a better  understanding  of  the 
SGEMP  response  of  complex  structures,  to  determine  how  well  the  SGEMP 
response  of  complex  structu'ces  are  represented  by  simple  bodies,  and  to 
determine  'now  much  structural  detail  must  be  included  in  the  model  to 
represent  well  defined  but  quite  complex  geometries.  Tne  purpose  of  this 
paper  is  to  report  ongoing  investigations  in  this  area,  and  in  paiticular 
to  address  (1)  representation  cf  complex  bodies  by  simple  spheres  and 
cylinders,  (2)  separability  of  the  internal  response  from  the  external 
response,  (3)  renresentation  of  segmented  bodies  by  continuous  bodies,  and 
(4)  response  ch;--,acceristi  cs  of  dumbbell  objects  ccnsisting  of  two  large 
bodies  separated  by  a boom. 

C-2.  COMPUTATIONAL  TECHNIQUE 

Ihe  ABORC  (arbitrary  body-of-revolution)  code  (Nef.  4)  used  in  this 
study  solves  the  complete  set  of  Maxwell's  equations  with  self-consistent 
electron  motion  for  axisymmetric  geometries.  Direct  finite-differencing 
of  the  field  equations  is  done  employing  generalized  coordinates,  and 
finite  "particles"  of  charge  are  followed  through  the  spatial  mesh  of 
zones  to  obtuin  currents.  Emission  of  arbitrary  energy,  angular,  soatial, 
and  time  distributions  of  currents  can  be  specified.  Randomizing  techniques 
for  the  en.ission  are  employed  for  efficient  representation  of  emission 
distributions.  Finite  conductivities  can  be  specified  representing  imperfect 
conductors,  and  dielectric,  structures  may  be  treated  by  specifying  proper 
dielectric  constants. 

T-E.  F.  Wenaas,  .S.  H.  Rogers,  and  A.  J.  Woods,  "Sensitivity  of  SGEMP 
response  to  Input  Parameters,"  IEEE  Conference  on  Nuclear  and  Space  Radiation 
Effects,  July  14-17,1975. 

2t.  N.  Delmer  et  al.,  "SGEMP  Phenomenology  and  Computer  Code  Development, 
DNA  3653F,  November  11,  197^. 

^E.  P.  Wenaas  et  al . , "Topics  in  SGEf4P  Analysis,"  IRT  document  INTEL-RT 
0001-080,  February  25,  1974. 

"^T.  A.  Tumolillo  et  al . , "Sk;met  Program:  Current-Injection  Predictions," 
3 volumes,  IRT  document  INTEL-RT  8121-007,  February  1975. 


The  free-space  boundary  condition  is  treated  by  enclosing  the  entire 
problem  in  a second  outer  conducting  body  of  revolution  rather  than  by 
utilizing  a radiative  boundary  condition.  Free-space  solutions  can  be 
obtained  by  moving  the  outer  boundary  out  so  the  clear  time  (the  time 
at  which  reflections  from  the  outer  wall  return  to  the  structure)  i.. 
larger  than  the  problem  time  of  interest. 

C-3.  MODELING  OF  EXCITATION  PARAMETERS 

All  geometries  have  been  placed  inside  a cylind-'ical  outer  enclosure 
of  30  m length  and  diameter.  This  provides  for  a clear  time  of  about 
100  nsec,  Thiis,  free-space  conditions  for  the  satellite  models  are 
simulated  for  a period  of  time  long  compared  to  the  assumed  incident 
photon  pulse,  modeled  as  o symmetric  triangle  with  a 10-nsec  rise  and 
fall  and  a 10-nsec  full  width  at  half  maximum  (FWHM) . The  lO-nsec  pulse 
used  in  these  studies  was  chosen  because  it  is  approximately  equal  to  the 
transit  time  for  light  around  the  skynet  satellite,  thus  allowing  for 
the  possibility  of  exciting  resonances  on  the  structures.  Results  obtained 
here  may  be  scaled  to  other  object  dimensions  and  pulse  lengths  usiiig 
appropriate  scaling  laws  (Ref.  2).  All  geometry  models  are  rotationally 
svnmetric , and  aK  surfaces  are  perfect  conductors. 

Electi'on  emission  curients  are  specified  from  the  various  surfaces 
of  the  bodies  representing  photo-electron  emission.  The  emitted  electron 
energy  spectrum  employed  ir.  ail  cases  is  shown  in  Figure  C-1.  The  angular 
distribution  of  emitted  electrcns  is  assumed  to  be  proportional  to  cos  6 
in  all  case.'"  '.'niform  spatial  distributions  of  emission  currents  are  assumed 
ir.  all  casts  except  the  Skynet  modet,  although  measurements  inaicate  a weak 
dependence  on  the  yield  weith  the  angle  of  incidence  of  the  photons  relative 
to  the  surface  (R.ji.  5)  . 

5m.  J.  Bernstein  and  X W,  Paschen,  "Forward  and  Backward  Photoemi ssion 
Yields  from  Metals  at  Various  X-ray  Angies  of  Inciden re , IEEE  Trans.  Nucl . 
Sci.  NS-20  (19731 . ' 
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Figursr  t.’-j  . The  electron  energy  snectrun'.  used  in  all  photo-en’issions 

here  is  the  revei-se  emission  from  A as  computed  by  the  POEM 
code.  'Ihe  photon  spectrum  is  a 50-keV  bi  cmsstrahlung  spectrum 
with  a 5-niil  Mylar  filter;  the  pr^^ponderance  of  photons  lie 
in  the  energy  range  from  7 to  50  keV  (Refs.  6,7). 

All  emissions  of  currents  are  '.Ime-phased  according  to  the  time  required 
for  the  photons  to  reach  the  emission  point.  The  calculations  have  been 
performed  for  two  fluence  conditions  wnich,  for  the  spectrum  indicated, 
result  in  peak  emission  currents  of  1.5  and  1500  amp/m  . Tne  lower  fluence 
corresponds  to  levels  at  which  the  fields  are  small  enough  that  no  significant 
electron  trajectory  perturbaticr.i  occur  and  the  problem  is  essentially 
linear,  while  high  fluence  indicates  significant  modification  of  electron 
trajectories  by  the  fields.  Results  presented  here  are  primarily  for  the 
high-fluence  regime  where  responses  are  nonlinear  and  effective  excitations 
are  shorter,  causing  higher -frequency  response. 

N.  Bradford,  "X-Rr.y- Inducer  Electron  Emission  II,"  lEEF  Trans 
Nucl  ■ Sci  ■ NS  - 20  1 1 9 '^3 J . 

■0.  N.  Bradford,  IEEE  Trans.  Nuc''. . Sci.  NS  - 19  (1972). 
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All  data  for  v/hich  time  histories  are  presented  are  averaged  over 
1 nsec  unless  otherwise  noted.  This  step  is  helpful  in  some  cases  for 
interpreting  results  due  to  numerical  hash  which  is  caused  by  the  particle 
nature  of  the  code. 

The  validity  of  some  of  the  assumptions  made  here  for  photo-emission 
specifications  with  respect  to  practical  considerations  can  ceitainly  be 
questioned  To  assutie  the  same  spectrum  for  forward-emitted  electrons 
as  for  ‘'a.ckward-emitted  elect^'ons  may  not  be  unrealistic  for  surfaces  such 
as  Mylar  thermal  blankets,  but  it  is  certainly  nor.  realistic  for  thick 
surfaces  or  where  significantly  different  materials  are  present.  To  use 
the  same  peak  intonsity  and  angular  distribution  from  the  sides  of  a 
cylinder  as  from  the  top  could  also  be  questioned.  However,  as  stated 
before,  the  purpose  of  this  paper  is  to  investigate  geometry  effects,  and 
we  have  endeavored  to  Keep  the  studies  a7id  results  uncluttered  with  effects 
and  numerou.s  carameter  variations  which  may  cloud  the  geometry  effects 
themsej.ves . Sensitivities  of  these  excitation  paramevers  are  treated  in 
the  companion  oaper  (Ref.  1) . 

C-4.  RbSPOHSE  CHAHACTERIZaTICN 

Wavs  tT  chaiacterite  the  response  of  an  object  iticlede  electric  and 
magnetic  fields,  surface  currents-,  charge  densities,  potentials,  etc. 
inasmuch  as  we  are  interested  in  the  response  of  the  structure,  tlie  fields 
at  the  surface  of  the  .structure  are  particularly  u.seful,  and  in  fact,  the 
normal  electiuc  field  and  tangertial  magnetic  field  at  the  surface  of  a 
perfectly  conducting  bedy  are  sufficient  to  specify  charge  dens:;  lies  and 
surface  currents.  In  this  })aper,  the  magnetic  fields  a.-e  used  pri.narily 
to  characterize  the  response,  although  it  must  be  remembered  that  the 
electric  fields  ere  just  as  important  in  many  cases. 

C-.S.  REPRIiSENfAi  IONS  OF  COM.h^.EX  EGDiHS  BY  srifLE  BODT  fS 

Tile  first  js.sue  to  bo  addrcs.sed  ir.  how  \;ell  complicated  satellite 
structures  can  be  represented  by  rather  simple  bodies 


revolution  such 


as  a cylinder  or  sphere.  This  question  is  particularly  relevant  because 
much  of  the  previous  modeling  of  the  SGEMP  response  of  structures  has 
been  performed  vith  these  simple  bodies  (Ref.  2) . For  purposes  of  this 
study,  we  compare  the  response  of  the  Skynet  satellite  (Ref.  4)  of  that 
of  a cylinder  and  sphere.  A simplified  representation  of  the  satellite 
is  shown  in  Figure  C-2,  along  with  a simple  cylinder  of  the  same  basic 
dimensions  and  a sphere  cf  radius  0.72  m having  a surface  area  equal  to 
that  of  the  cylinder.  The  regions  of  emission  on  all  three  structures 
are  represented  by  the  shaded  areas  on  the  surfaces.  Higher  electron 
emission  from  regions  of  the  satellite  plated  with  gold,  are  indicated  by 
lOX,  while  regions  of  the  satellite  with  lower  emission  are  indicated 
by  0.2X.  No  emission  is  assumed  from  the  inside  of  the  outer  wall  because 
this  surface  is  coated  with  low-Z  material  which  emits  fewer  electrons. 

The  outer  wall  is  modeled  as  an  electrically  isolated  surface  to  represent 
its  being  held  in  place  by  dielectric  braces.  In  reality,  the  outer 
structure  is  electrically  connected  by  means  of  solar  cell  cables,  but 
these  cables  cannot  be  modeled  by  a body  of  revolution  without  disastrous 
results  for  the  structural  response.  The  problem  of  two  bodies  connected 
by  a cable  or  thin  rod  is  treated  in  subsequent  sections. 

The  regions  of  emission  on  the  outer  surface  for  the  three  bodies  aie 
chosen  so  that  the  emission  surface  areas  are  approximately  equal.  Half- 
sphere emission  is  chosen  for  the  sphere  to  generate  worst-case  surface 
curients  for  the  spherical  model.  This  constraint,  along  with  that  of 
eequal  emission  surface  areas,  requires  an  unrealistic  emission  pattern 
from  the  cylindrical  object  in  which  emission  occurs  over  the  top  half  of 
the  structure  only.  Thus,  it  is  evident  that  the  sphere  has  at  least  one 
shortcoming  in  modeling  a cylindrical  structure. 

Results  ill  terms  of  the  magnetic  fields  at  the  center  of  the  side 
surfaces  of  the  three  models  for  the  high-fluence  case  are  shown  in  Figure 
C-.'i.  'Hic  prcdi -led  )H-ak  currents  for  tlie  three  structures  are  similar, 
although,  the  simple  iiio'.lels  f a ' 1 to  reproduce  the  resonant  beliavior  exhib- 
ited bv  tin-  more  eoiii])]  ; cated  body.  In  general,  the  modeling  of  the  exter- 
ior suri'acc  ciirren’.s  !iv  simple  !;odies  appears  to  be  reasonable. 
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EMISSION  REGIONS 


Geometry  models  investigated  with  ABORC  for  determination  of  validity  of 
representing  satellite  structures  with  simple  cylindrical  and  spherical  objects 


Figi’re  C-3.  Comparison  of  magnetic  fields  at  the  sides  of  a sphere,  cylinder, 
and  Skynet  model.  The  peak  emission  current  is  1500  amp/m2, 
corresponding  to  a highly  space-charge-limited  solution. 


'fhe  response  of  the  complex  body  is  shown  in  Figure  C- 1 iit  a second 
point  within  the  reentrant  portion  of  the  structure  which  has  no  parallel 
on  the  simple  bodies.  The  response  at  this  point  is  quite  different  from 
the  responses  exhibited  by  the  simple  bodies.  Thus,  we  conclude  that  simple 
bodies  can  reasonably  model  the  external  response  at  least  in  terms  of  peak 
currents,  but  that  they  cannot  model  features  introduced  by  complex  struc- 
tures such  as  reentrant  bodies. 
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Figure  C-4.  Magnetic  field  between  equipan.  box  and  outer  wall  for  Skynet 
satellite  model.  The  ak  emis:.inn  current  i 1500  amp/m^, 
corresponding  to  a higi  S|'  e-charge-limited  solution. 

C-6.  SEPARAJ’TLI'/  Qi-  TNT.fRNAL  ANo  XTi  ^ \L  RESPONS) 

TT.e  modeling  of  the  fl.  ^ response  ot  i,  lex  bodies  us  ally 

performed  by  separating  nhe  n onse  of  the  i.nt,;rn,  p ortion  of  he  problem 
from  that  of  the  external  oortinn.  Thus,  th  details  oi'  the  internal 
nonlinear  response  are  assumed  i be  unpertuiucd  by  external  emission  and 
vice  versa.  Leakage  currents  between  the  nterior  and  exterior  are  con 
sidered  although  they  are  assumed  to  repress  a sinall  perturbation,  and 
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the  two  solutions  are  assumed  to  be  additive.  Obviously,  this  assumption 
is  valid  in  the  low-level  linear  regime  where  field  coupling,  even  if 
significant,  will  not  affect  the  electron  motion.  Thus  the  separability 
issue  is  treated  here  for  the  high-fluence,  nonlinear  response  only. 

The  geometry  of  Figure  C-S  has  been  chosen  because  such  a configura- 
tion, with  the  relatively  large  9-cin  gap,  v.as  thougl.t  to  maximize  the 
coupling  between  tiie  external  and  internal  response,  thus  tending  to  rep- 
resent a worst  case.  In  fact,  it  will  be  shown  in  the  following  section 
that  the  coupling  of  the  two  responses  is  much  less  severe  for  a similar 
object  without  the  center  rod.  A gap  of  f cm  is  larger  than  most  gaps  in 
actual  s.itellites,  <aiid  this  again  tends  to  be  worst-case. 

The  electric  and  magnetic  fields  were  computed  at  the  f>ui  points 
shown  in  Figure  C-5  for  emission  from  the  surfaces  irdicared  by  the  .shaded 
and  dashed  lines. 

Three  cases  were  considered,  including  (1)  exte.’'nal  -emission  only, 

(2)  internal  emission  only,  and  (.3)  simultaneous  internal  and  external 
emission.  The  resulting  exterior  magnetic  field  at  point  1 i?  shown  in 
Fig’jre  C-6  for  exterior  emission  only  and  interior  emission  only.  The 
resulting  interior  magnetic,  field  at  point  2 is  shown  in  Figure  C-7,  again 
for  exterior  emission  only  and  interior  emission  only.  It  i.s  evident  from 
these  results  that  emission  from  the  external  surfaces  causes  currents  to 
flow  on  the  interior  which  are  on  the  same  order  as  those  produced  by 
internal  emission,  but  that  the  converse  is  not  true.  'Fliat  is,  currents 
produced  by  emission  on  the  in.side  do  not  cause  large  structural  currents 
on  the  outside. 

Tne  magnetic  fields  are  n'-t  the  entire  story,  however.  Perturbations 
in  the  nonlinear  response  arise  from  perturbations  in  the  electric  field 
rather  than  the  ma  ir.etic  field  (e.g.,  the  nonlinearity  in  electron  trajec- 
tories is  caused  primarily  by  the  electric  field,  not  the  magneti:  field) . 
Thus,  the  strong  coupling  of  the  outside  current  to  the  inside  does  not 
necessarily  mean  that  nonlinear  internal  response  cannot  be  computed 
separately  from  the  external  response.  The  peak  fields  are  summarized  in 
Table  1,  and  the  fields  at  two  points  in  time  in  Table  2,  where  it  is 
evident  that  the  electric  field  coupling  from  the  inside  to  outside  is 
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Table  C-1 


PHAK  FIHLDS  AT  POINTS  INDTCATHD  IN  FIGURF  C-5  FOR  (1)  INSIDE  EMISSION  ONLY, 
(2)  OUTSIDE  EMISSION  ONLY,  (3)  SUMMATION  OF  INSIDE  ALONE  AND  OUTSIDE  ALONE, 
AND  (4)  SIMULTANEOUS  EMISSION  (UNITS  ARE  amp/m  OR  V/m) 


Inside  Outside  Superposil ion  ot  Inside- 
Emission  Emission  Only  and  Outside-Only  Simultaneous 


Field  Points 

Only 

Only 

"l 

3.1 

15 

-104 

-80 

"3 

-116 

-75 

-3,7 

-3.2 

4 

4 

5 

^1 

2.0x10 

1 . 1x10 

5 

4 

E„ 

-I . 1x10 

-2.0x10 

JL 

4 

4 

2.0x10 

2.0x10 

-2.2xl0'^ 

4 

-3.1x10 

Solutions  Emission 


18 

20 

-184 

-128 

-191 

-136 

-6.7 

-4.2 

1.3x10^ 

1 . 1x10 

-1 . 3x10^ 

■1 . 1x10 

4.  Oxio'^ 

2.6x10 

-5.3x10"'’ 

-3.4x10' 

Table  C-2 

FIELD  VALUES  AT  POSITIONS  INDICATED  FOR  (1)  INSIDE  EMISSION  ONLY, 
i>)  OUTSIDE  EMISSION  ONLY,  (3)  SUMMA'ilON  OF  INSIDE  AI.ONE  AND  OUTSIDF  ALONE, 
AND  (4)  SIMULTANEOUS  EMISSION;  FIELDS  IND1CA''ED  AT  10  nscc, 

TUL  APPROXIMATE  TIME  OF  PEAK  OF  INTERNAL  RESPONSE  (UNITS  ARE  V/m) 


F^eld  Points 

Inside 
Emiss ion 
Only 

Outside 

Eniissior. 

Only 

Superposition  of  Ir.side- 
Only  and  Outside-Onl/ 
Solutions 

Simultaneous 

Emission 

Time  - 10  nsec 

E^  (outside) 

E„  (inside) 

1 . 1x10''’ 

.S 

1 1x10 

1 . 1x10''’ 

-1.1x10^ 

-7.7x10'"’ 

-1.2x10''’ 

-!  . Ixio'’ 

(insiJe) 

1.4x10'’ 

-7.7x10'^ 

X 

o 

6 . 5x10'’ 

3 

3 

3 

E^  (outside) 

2.6x10 

-6. 9x10 

■ 4.3x10 

r- 

-4.5x10 

rime  = 30  nsec 

E^  (outside) 

E,,  (inside) 

5/7x10'^ 

4 

2.9x10 

3.5x10)'’ 

2.7x10“’ 

■7..‘.xUr'^ 

-2.3x1  O'"’ 

-9. 8.x  10'^ 

-5  , 3x10'’ 

1 

E..  (’nsidej 

4 

1 ..‘.xlO 

4 

1 . 7x10 

4 

3.  dxlO 

i.  9x1  o'’ 

E‘  (ofi'LS  Lde  ) 

- 1 . 2XiO''' 

4 

-2.4x10 

- .i  . (iV  1 0 

•2. 4x1  o’ 

tvl 


weak,  and  vice  versa.  Thus,  there  is  every  reason  to  believe  that  the  two 
nonlinear  solutions  can  be  computed  separately.  This  concluc.ion  is  demon- 
strated in  Figures  C-8  and  C-9,  where  the  solutions  obtained  by  adding 
the  response  due  to  external  emission  alone  to  the  response  from  internal 
emission  al-'ne  are  compared  to  the  response  with  simultaneous  emission 
from  the  inside  and  outside.  Thus,  the  internal  and  external  responses 
appear  to  be  separable  in  that  the  notilinear  responses  may  be  computed 
independently.  However,  the  external  response  is  in  general  not  isolated 
from  the  internal  response,  and  the  field  leakage  must  be  considered. 


1 « 


Figure  C-S.  Geometry  utilized  in  separability  study.  Four  field  points 

arc  considered  a.s  indicated.  The  radial  position  for  H-fieids 
is  0,34  m and  for  £• fields  0.41  m. 
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Figure  C-6.  Magnetic  field  at  point  1 outside  the  cylinder  obtained 
with  exterior  smis  ion  alone  compared  with  the  response 
at  the  same  point  produced  by  interior  emiss'on  alone 
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Figure  C-7.  Magnetic  field  at  point  2 insvde  cylinder  obtained  with 
internal  emission  alone  compared  with  response  at  same 
point  produced  by  external  emission  alone 
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Figure  C-8.  Magnetic  field  at  point  1 outside  the  cylinder  obtained  by  add- 
ing the  response  from  inside  emission  alcne  to  the  response  from 
outside  emission  alone  compared  to  the  response  at  the  same  point 
obtained  by  simultaneous  emission  from  the  inside  and  outside 
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Figure  C-R.  Magnetic  field  at  point  2 inside  the  ivlinder  obtained  by  adding 
the  I'Csponse  from  inside  emission  alcne  to  che  response  from  out- 
side emission  alone  compared  to  the  response  at  the  same  point 
obtained  by  simultaneous  emissiori  from  tlie  inside  and  outside. 
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C-7.  SEGMENTI-D-BODY  EFFECTS 


Sarellite  bodies  have  generally  been  represented  by  continuous  bodies 
of  revolution  for  purposes  of  SGEMP  calculations,  although  many  of  the 
spacecraft  manufacturers  have  designed  these  systems  with  segmented  panels, 
as  shown  in  Figure  C-lOa.  Although  this  type  of  segmentation  cannot  be 
modeled  wi tii  codes  limited  to  rotational  symmetry,  we  can  gain  insight  into 
the  response  of  segmented  bodies  by  considering  segmentation  in  the  axial 
direction,  as  shown  in  Figure  C-lOb,  ri.o  figure  without  a connecting  rod 
simulates  truly  isolated  panels,  while  the  configuration  witli  a connecting 
rod  simulates  an  internal  connection  between  panels  such  as  a solar  cell 
cable . 

Computations  were  performed  with  the  three  bodies  shown  in  Figure 
C-iOb  with  equal  emission  from  the  top  half  of  each  cylinder.  Peak  field 
values  are  shown  in  Table  C-3  for  the  three  geometries.  The  table  reveals 
very  small  differences  bet\-een  outside  peak  field  values  due  to  segmenting 
or  interior  return  paths.  The  slight  difference  in  the  peak  electric  field 
at  the  back  in  the  separated  cylinder  case  is  due  to  the  lack  of  charge 
transfer  to  the  bottom  half  of  the  cylinder  in  That  case.  A very  signifi- 
cant effect  on  the  interior  magnetic  field  is  seen,  due  to  large  curr'cnts 
on  the  interior  of  the  body  flowing  through  the  rod. 

The  external  magnetic  field  near  the  top  of  the  cylinder  is  shown  in 
Figure  C-11  for  the  case  of  high  space-charge- 1 imi T ing . It  is  quite  evi- 

dent that  the  response  is  essentially  the  same  for  these  objeets.  Similar 
results  also  hold  for  the  magnetic  fields  at  other  points  on  ttie  external 
surfaces . 

Differences  in  electric  fields  on  the  lower  half  of  the  cylinder  should 
be  evident  in  the  isolated  cylinder  as  compared  ro  the  conducting  cylinder 
due  to  the  difference  in  charge  transfer.  (The  late-time  electric  field  on 
the  lower  naif  of  the  solid  bodv  should  be  larger  than  on  the  isolated  body 
because  t r.ir.s ferred  charge  can  be  redistributed  on  the  solid  body  while  it 
cannoi  redistribute  in  the  isolated  case  except  by  capacitive  coupling,) 
Ixpected  differences  in  the  late-time  electric  fields  at  the  bottom  of  the 
cylinder  are  evidiut  in  Figure  C-12. 
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a.  TYPICAL  SATELLITE  SEGMENTATION 
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b.  SEGMENTED  BODIES  UTILIZED  FOR  EXAMINING  SEGMENTED  BODY  EFFECTS. 
S'-IPLC  CYLINDER  AT  RIGHT  IS  COMPARED  WITH  RESULTS. 

RT-11863 


Figure  C-lU(a),(b).  Typical  and  |.>resenv  1 v calculable  segmented 

bodies  in  nonlinear  SGl.MP 


Table  C-3 


PKAK  FIF.LD  VALUES  AT  F^OSiTIONS  INDICATED'’* 

ON  A SIMPLE  CYLINDER  AND  SEGMENTED  CYLINDERS  OBTAINED 
AT  HIGH  FiAJENCE  (UNITS  ARE  amp/m  OR  V/m) 


Field 

Solid 

Cylinder 

Hoi  low 
Cylinder 
Segmented 

Segmented 
Hollow  Cylinder 
Connected  by  Rod 

13 

14 

14 

= '2 

0 

5 

-SO 

0 

6 

-75 

^’4 

-4.9 

-3.2 

-3.2 

4 

5 

5 

3 

^1 

1 .1x30 

1.1x10 

d 

1 . 1x10 

4 

1.2 

0 

-2.8x10 ■ 

4 

-1.8x30 

4 

0 

-1  . 2x10 

2.0x10 

4 

4 

4 

-3.1x10 

-1  .SxlO 

-3. 1x30 

^Point:;  2 a-'d  3 arc  inside  the  top  and  bottom  walls  of 
the  cylinders,  respectively.  Radial  position  for  il 
values  is  0.34  m,  for  E values  0.41  m. 
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MAGNETTC  FIELD  (amp/m) 


Figure  C-11.  Exterior  magentic  fields  at  top  surface  (radius  0.34  m) 
versus  time  for  simple  and  segmented  cylinders  shown  in 
Figure  C-IO 


RT-n864  TIME  (nsec) 

Figure  C-12.  Exterior  electric  field  at  point  4 (radius  0.41  m)  outside 
the  cylinder  at  bottom  for  simple  and  segmented  cylinders. 
Peak  emission  current  of  1500  amp/ni2,  corresponding  to  high 
space-charge-limiting,  was  used. 


internal  response  of  the  isolated  r-  Under  is  rr,arkedly  different 
fi'om  the  response  of  the  cylinder  with  the  rod,  however.  This  is  not  sur- 
prising in  as  much  as  current  is  driven  internally  because  of  the  transfer 
of  charge  between  segments.  Large  low-frequency  oscillations  are  evident, 
representative  of  an  LC  circuit  consisting  of  the  inductance  of  the  rod 
and  capacity  of  the  segments. 

Ir  summary,  the  modeling  of  the  gross  external  currents  of  a segmented 
body  can  be  well  represented  by  a body  of  revolution,  but  in  general,  the 
response  within  tlie  object  is  lest  by  such  simple  modeling. 

C-S.  DUMPBELL  GEOMETRY 

Another  facet  of  modeling  which  has  been  exp.lored  is  a geometrical 
configuration  consisting  of  two  relatively  large  bodies  separated  by  a 
thin  red  or  boom.  Such  a configuration  might  correspond  to  the  case  of 
solar  paddles  extending  fiom  a satelrite  body. 

The  dumbbell  geometry  shown  in  Figure  C-13  was  chosen  for  this  study. 
The  dimensions  of  the  mam  bray  are  the  same,  as  those  of  the  isolated 
cylinders  treated  previously.  Uodit.  lengths  of  1 and  2 m were  considered, 
with  a smaller  body  attached  to  the  end.  The  boom  ladius  was  held  constant. 
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Figure  C-  13,  Dumbbell  geometry  calculated  v/ith  ABORC  for  study  of  sfcprrate 
bodies  connected  by  booms.  Boom  lengths  of  1 and  2 m are 
considered . 
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Calculations  were  performed  for  emission  from  the  top  and  side  of  the 
siTiall  cylinder  and  from  the  top  half  of  the  large  cylinder,  as  shown  in 
the  figure.  Thus,  the  emission  from  the  larger  cylinder  is  identical  to 
that  used  in  the  isolated-cylinder  studies. 

Figure  C-  14  shoi/s  the  magnetic  field  responses  at  the  midpoint  of  the 
larger  cylinder  of  the  dumbbell  compared  to  results  obtained  by  treating 
the  larger  cylinder  separately.  The  responses  at  this  point  are  strikingly 
similar,  indicating  that  the  boom  and  small  attached  object  do  not  signi- 
ficantly perturb  these  fields  on  the  side  surface  of  the  large  cylinder  at 
this  high-fluence  level. 

The  response  of  the  boom  is  shov/n  in  Figure  C-15  for  boom  lengths  of 
1 and  2 m.  It  is  evident  that  the  response  on  the  boom  is  significantly 
different  from  tbe  response  on  the  body.  The  large  resonances  appear  to 
be  characteristic  of  ao  LC  circuit  consisting  of  the  inductance  of  the 
dumbbell  and  the  capacity  of  the  two  bodies.  The  resonant  frequency 
decreases  with  increasing  boom  length,  as  expected. 

Peak  il-fields  in  the  boom  ^re  significantly  higher  than  on  the  sides 
of  the  body,  but  the  total  peak  currents  flowing  on  the  boom  are  less  than 
on  the  large  body: 

Iboom  “ ^ ^TTr  - l.>5  x 2Tr(0.07}  = 57  amps  , 

^boom  ~ ^ 4;t(0.7)  ~ i80  amps  . 

Thus,  it  is  evident  why  the  presence  of  the  boom  does  not  significantly 
perturb  rhe  res])on.jc  of  the  larger  cylinder.  The  currents  on  the  boom 
are  lower  than  on  the  surface  of  the  body  due  to  effective  cancellation 
from  currents  emitted  to  infinity  from  each  body. 

At  the  higher  fluence  levels  treated  here,  the  sj)ace-charge- 1 imi  ted 
responses  of  each  body  appear  to  be  somewhat  decoupled.  As  indicated  in 
Table  C.-4  , the  peak  magnetic  fields  do  not  change  s i gn  i t i can  t ly  as  the  two 
bodies  are  moved  further  apart,  indicating  that  the  trajectories  of  elec- 
trons emitted  from  the  two  bod  es  are  not  altered  by  the  presence  of  each 
other.  In  the  case  of  lower  siiace-chargc- 1 imi t ing , where  a large  number 
of  electron  trajectories  would  reach  oat  to  distances  on  the  order  of  the 
separation  'f  bodies,  however,  we  would  expect  much  more  coui)ling. 


Table  C-4 


SEGMENTED- BODY  EFFECT;  PEAK  N5A.GNETIC  FIELD  VALUES 
AT  POINTS  INDICATED  AND  OSCILLATION  PERIODS  FOR  A CYLINDER 
COMPARED  TO  A DUMBBELL  CONFIGURATION  AT  HIGH  FLUENCE 


C-9.  SUMMARY  AND  CONCLUSIONS 


The  sel  f-con;;  i stent , fully  dynamic  ABORC  code  for  rotational  ly  sym- 
metric bodies  has  been  employed  on  a number  of  bodies  to  investigate  the 
effects  of  approximations  of  geometries  in  analyzing  the  SGEMP  response 
of  a system.  We  arrive  at  the  following  tentative  conclusions  based  on 
this  relatively  small  number  of  geometrical  studies. 

1.  The  response  of  simple  bodies  such  as  cylinders  and  sphere:^ 
is  similar  tc  more  complex  geometries  in  terms  of  peak  cur- 
rents and  response  times. 

2.  The  inside  and  outside  responses  of  an  object  can  be  sep- 
arated to  first  order  in  many  cases,  even  though  both  solu- 
tions by  themselves  are  highly  nonlinear  and  the  leakage 
current's  arc  relatively  large.  'ior.ie  interaction  of  the 


I internal  and  external  responses  was  seen*  and  that  is 

I undoubtedly  geometry -dependent . 

{ 3.  The  external  response  of  s highly  segmented  conducting  body 

‘ is  similar  to  the  response  of  a smoothly  connected  bedy  of 

!;  revolution.  The  internal  response  can  be  considerably  dif- 

' ferent  when  a conducting  path  exists  between  the  segments, 

i The  response  then  contains  a long  low-frequency  oscillation 

I similar  to  the  response  of  an  LC  circuit  consisting  of  the 

I inductance  of  the  interconnecting  rod  and  the  capacity  of 

i the  segments. 

I 4.  The  space-charge-iimiting  characteristic.s  of  two  large 

objects  separated  by  a distance  large  compared  to  the  dimen- 
sion  of  the  space-charge  barrier  are  relatively  unperturbed 
? by  the  presence  of  each  othe^. 


D-1.  INTRODUCTION 


A question  which  must  be  considered  with  every  calculation  made  by 
a finite-diffsrence/particle-mover  code  is  how  sensitive  the  results  are 
to  the  numerical  grids  and  particle  statistics  employed  in  the  problem. 

Th?s  issue  is  particularly  important  when  calculations  are  being  compared 
where  only  minor  changes  in  input  are  being  considered  to  determine  effects 
on  response.  In  this  appendix,  some  results  of  various  grid  sensitivity 
investigations  carried  out  during  the  course  of  several  projects  are 
presented.  Tiie  goal  is  to  demonstrate  how  to  evaluate  ABORC  calculations 
for  numerical  grid  sensitivity.  Some  example  calculations  are  shown,  and 
parameters  useful  in  measuring  calculational  quality  are  discussed. 

This  discussion  is  an  outgrowth  of  several  project'-,  and  therefore, 
the  examples  may  be  mildly  disjointed.  The  principles  discussed  apply 
to  a broad  spectrum  of  conditions,  however. 

D-2.  ABORC  GRIDS 

Grids  which  must  be  sjiecified  in  ABORC  calculations  are; 

Spatial  coning. 

Time  steps. 

Energy  bins. 

Angular  distributions. 

Emission  spatial  zones. 

The  first  three  grid  types  are  considered  here.  Results  of  variations 
of  the  la.st  two  input.s  on  SGEMP-iEMP  calculations  can  be  found  in  Refer- 
ences 1 and  2. 

^E.  P.  Wenaas,  S,  Rogers,  and  A.  ,J . Woods,  "Sensitivity  of  .SGP.MP 
Response  tc  Input  Parameters,"  TEI'E  Trans.  Nuc . Sci.,  Vol  . NS- 22 , Dec.  1975, 
p.  2362. 

'E.  f.  dePlomb  and  A.  J.  Woods,  "TEDiEM-RZ  and  RG;  Two- D imen;.  i ona  i 
Tine- Dependent  lEMP  Computer  Codes,"  DNA  3140F,  Mar.  10,  1973. 


Considerations  which  must  be  made  in  choosing  ABORC  grids  are  summar- 
ized below. 

Resolution  of  object 
Resolution  of  time  response 
Pulse  shape 
Electron  velocity 
Clear  time 
Field  gradients 
Statistical  noise 
Computer  run  time 
Computer  memory 

Object  resolution  can  require  fine  grids  if  effects  of  re-entrant 
bodies  are  being  considered,  or  if  body  shapes  are  unnatural  for  the 
cylindrical  coordinates  employed  in  ABORC,  or  if  pulse  lengths  are  short 
compared  to  object  dimensions.  For  simple  cylinders,  however,  body 
resolution  is  usually  not  a restriction.  Time  response  resolution 
requirements  also  must  be  considered,  particularly  if  the  problem  is 
very  dynamic.  At  least  5 or  6 grid  points  are  required  for  accurate 
transmission  of  wavelengths.  Shorter  wavelengths  are  distorted  or  do 
not  propogate  through  the  gird  at  all  (Ref.  3).  Obviously,  the  pulse 
shape  must  be  considered  in  choosing  time  steps  and  spatial  grids  if 
it  is  short  or  has  fine  structure  which  should  be  resolved.  Electron 
velocity  is  important  in  that  the  distance  travelled  by  tliC  electron; 
in  a time  step  should  be  roughly  compatible  wirii  grid  spacing  and  field 
gradient  distances.  Clear-time  requirements,  or  time  for  radiation  to 
be  reflected  back  to  the  inner  object  from  the  outer  tank  walls,  generally 
determine  the  outer  cylinder  size  because  ABORC  does  not  have  a free- 
space  boundary  condition.  A large  clear  time  results  in  large  outer 
zone  and  charge  zone  size  changes,  especially  if  small  object  detail  or 
space-charge  barriers  are  being  ie.-,olved. 

Distances  over  which  fields  change  significantly  shoulii  also  be 
considered  in  choosing  space  and  time  zones.  ^ny  gross  discrepancy  in 

Boris,  "Relativistic  Plasma  Simul  at  i on- Opt  i mi  zat  ion  of  a 
Hybrid  Code,"  Proc . fourth  Conference  on  i'.uinerical  Simulation  of 
Plasmas,  NRL,  IQ^O. 
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which  resolution  of  a space-charge  barrier  much  smaller  than  ;he  minimum 
spatial  zone  size  is  attempted,  for  example,  is  bound  to  result  in  grid- 
dependent  results.  Computer  memory  available  limits  the  code  to  100  by 
100  spatial  zones  at  present. 

Statistical  noise,  due  to  finite  particles  of  charge  crashing  through 
the  spatial  mesh,  decreases  with  the  number  of  particles.  Computer  run 
time  usually  limits  the  number  of  particles,  however.  The  time  required 
for  both  particles  and  fields  is  shown  below  in  7600  central  processor 
(CP]  seconds . 

-4 

Particles:  3 x 10  sec/particle  step 
Fields;  10  ^ sec/zone  step 

These  numbers  translate  into  about  5 minutes  of  CP  time  for  a calcu- 
lation with  an  average  of  3000  particles  followed  for  300  steps  through  a 
spatial  mesh  of  70  by  40  axial  and  radial  zones.  Tlie  piesent  cost  for 
such  a calculation  at  a government  installation  is  about  $100.  While  the 
field  solution  portion  of  the  code  is  programmed  very  efficiently,  very 
little  optimization  of  the  particle-pusher  coding  has  been  performed  to 
date.  Present  run  times  could  be  reduced  considerably  through  efficiency 
measures . 

Obviously,  the  large  number  of  considerations  in  choosing  grids 
prohibits  writing  down  an  exact  formula  for  grids  which  always  work 
withovit  resulting  In  unreasonable  computer  time  requirements.  Also,  to 
even  suggest  that  all  of  the  above  are  impoitant  in  producing  grid- 
independent  results  is  somewhat  speculative.  For  example,  the  requirement 
of  limiting  the  distances  over  which  particles  travel  in  each  time  step 
to  dimensions  small  or  comparable  to  distances  over  which  fields  acting 
on  them  change  considcr.ibly  seems  intuitively  very  reasonable.  The  })roblem 
is  so  complicated  in  a case  with  many  particles  ;it  different  energies  that 
results  might  smooth  out  atid  give  correct  responses  even  tlunigh  trajectories 
of  individual  jiarticles  might  be  grossly  different  from  reality.  'Pie  best 
way  to  know  for  cert;iin  is  to  vary  grids  while  holding  physical  properties 
of  a calculat  ion  const.ant  , and  compai-e  responses.  In  the  following  solutions, 
results  of  some  of  tliese  grid  v.ariations  are  shown,  and  discussion  of  some 
of  the  parameters  which  determine  numerical  quality  of  the  calculations  is 
g i veil . 
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Table  D-2 


NUMERICAL  GRID  SPECIFICATIONS  FOR 

MEDIUM- 

AND 

HIGH-FLUENCE 

ZONING  STUDIES 

Medium-Fluence 

High- 

Fluence 

A 

B 

A 

B 

Minimum  zone  size  (m) 

0.005 

0.02 

0.002 

0.005 

Number  of  emission 

zones 

3 

7 

4 

4 

Energy  bins 

Time  step  (sec) 

15 

15 

10 

10 

Particles  1 

.00x10'^*^ 

1.00x10' 

10 

2,5x10'^^ 

2.5x10'^^ 

Light  1 

. llxio'^  ^ 

3. 33x10' 

10 

- 12 

4.17x10 

8.53x10'^^ 

Number  cf  particles 
in  system  at  10  nsec 

2214 

2273 

1258 

1520 

Total  CPU  time  to 

10  nsec  (sec) 

M 

3() 

229 

160 

lixact  axial  zones  above  the  emission  surface  can  be  seen  in  Figures 
0-2  and  D-a  (similar  variations  in  radial  zones  were  undertaken  simultan- 
eously), along  with  effects  on  normal  electric  fields  in  tlio  space-charge 
r-ifrier  when  zone  sizes  are  varied.  Figure  l)-2  is  for  moderate  SCiL,  while 
i’igure  O-.l  is  for  very  high  SCI,.  Notice  the  much  steeper  gradients  and 
higiier  fields  in  the  liigh-SCL  case.  It  is  obvious  that  the  particles  will 
see  a considerably  higher  i.'lectric  field  right  near  ttie  emission  face  with 
the  fine  zoning  than  with  the  coarser  zoning.  Inis  is,  also  true  in  rlie 
medium-StiL  case,  only  mucti  less  pronounced.  It  would  be  reasonable,  there- 
fore, to  ex]H'ct  zone  size  sensitivity  in  the  liigh-SCI.  case  to  be  more  pro- 
nounced than  i I'  the'  iiiodei'ate- SCI,  case.  .Also  notice  that  much  coarser  zcnies 
can  be  employed  at  lower  fluences  witiiout  I'esultant  unreasonable  changes 
in  fiekls  from  giid  ju'int  to  g.rid  point. 

Time  histories  of  fields  al  sevei'al  positions  around  liie-  object  ai'e 
fouiKi  in  Figures  b-  1 aiul  It-.S,  Fig.ure  l)-  l is  tor  medium  i'l.tnice  a.  d D- S 
for  biiTi  fluoiK'e.  i Ik'  cuin(.'S  are  smoo.lied  with  a U.S-nsec  time  constant 

in  the  lim'd  i uiii- :1  uenci-  i.ase'  and  O.J  lu-ec  in  tlu'  h i gb  - f 1 ui'rr.'e  ci.se.  The 


D-3.  aeSULTS  OF  GRID  STUDIES 


Fluence  Dependence  of  Grid  Requirements 

Of  the  considerations  mentioned  above  for  choosing  grid  sizes  in 
ABORC  calculations,  the  resoluti^  n of  field  gradients  has  been  one  of 
the  most  consistent  problems  in  pi-oducing  grid-independent  responses. 
Gradients  are  not  generally  a problem  for  low-fluence  conditions  in 
SGEMF . 

As  fluence  is  increased  and  the  electrons  are  turned  back  by  fields, 
very  large  charge  densities  can  build  up  close  to  the  object's  emitting 
surfaces.  The  distance  over  which  the  electric  field  changes  by  a factor 
of  10  can  be  as  little  as  2 or  3 cm  for  fluences  and  spectra  of  interest, 
while  the  object  may  be  3 m long.  Such  li .ge  differences  in  zone 
requirements  to  describe  both  the  space-charge  barrier  and  the  object 
cause  concern  over  the  spatial  zone  dependences  of  results. 

To  ir-estigate  spatial  zone  sensitivities,  calculations  were  per- 
formed in  which  spatial  grids  were  changed  while  holding  all  other 
parameters  constant.  The  comparisons  were  made  at  intermediate  and  high 
SCL.  Physical  problem  co.iditions  are  listed  in  Table  D-1  and  numerical 
grids  in  Table  D-2.  Electron  emis.^ion  results  were  obtained  from  QUICKE2. 

Table  D-1 

TEST  PROBLEM  DESCRIPTIONS 
FOR  ZONING  STUDIES 


Problem  Geometry  Cylinder  with  length  - diameter  = 3 m 
Clear  Tine  Greater  than  simulation  time 

Elect rcii  Emission 

Emission  Surface  ; lop  ♦ half  side  urdform 
Spectrum:  see  Figure  D-1 

Angular  distribution:  cos  Q 

Time  history:  medium  fluence:  10-nsec  rise,  40-risec  fall 

2 

high  fluence;  sin  pulse  with  l/-nscc  rise 


medium- fluence  curves  show  almost  complete  indeoendence  of  the  grid  rhanj'e, 
whereas  substantial  differences  arc  seen  in  the  high-lev( 1 case,  especially 
for  field  quantities  evaluated  in  regions  where  many  parti  Ics  are  present. 
The  curves  permit  several  conclusions  to  be  drawn  I'egarding  spatial  zone 
requirements  in  ABORC  cal  eolations’ 

• Spatial  grid  requirements  are  fluence-dependent . 

• Time  histories  of  fields  near  the  object  are  virtually  spatial- 
zone-  independent  for  zoning  in  which  electric  fields  fall  off 
less  than  a factor  of  two  for  each  zone  (see  Figures  D-2  and 
D-4) . 

• Spatial  zone  dependence  of  results  is  observed  for  electric 
fields  falling  off  by  a factor  of  three  in  one  zone  (Figures 
D-3  and  0-5) . 

• Sensitivity  of  response  to  spatial  zoning  is  greatest  in  regions 
close  to  the  space-charge  barriers  (Figure  l)-5a  versus  Figure 
[)-5d)  - i.e.,  where  many  particles  of  charge  are  present. 

RKFKRF.NCHS 
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(V/m) 


DISTANCE  FROM  EMISSION  FACE  (m) 


D-2  !'i  sir  ihut  ions  of  normal  electric  field  near  peak  of  en.issiop 

current  pulse  for  different  S[)atial  zonings.  All  field  points 
near  the  emission  face  arc  plotted. 


(V/m) 


V AZ  = 0.005 


AZ  = 0.02  m 


qpc^<^^ 


h* — 0./5  m 
1 E ALONG  THIS  LINE 


- 8 X lo"^  amp/m^ 


TIME  - 10  nsec 
^RISE  " 


■TO  EMISSION  FACE 


SIGN 

CHANGE 


0.001 


RT-13666 


DISTANCE  FROM  EMISSION  FACE  (m) 


Figure  D-3.  Distributions  cf  normal  electric  field  near  peak  of  emission 

curren'  pulse  for  different  spatial  zonings.  All  field  points 
aie  plotted. 


I- i sure  |)-4r.  M;ii’,netic  field  lime  hi-. lories  at  side  of  cylinuci'  1 . (i  m 
from  ;<.!))  for  medium  SiiL,  obtained  with  nuninium  'qtatial 
zoniji).’,s  i.)f  O.Od  and  (,) . OOS  mi 


I ! 
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(V/m) 


t^T-13671  (nsec) 


I'igure  D .Sa.  [iltH'trir  field  time  histories  at  top  of  cylinder  at 
radius  0.7S  m for  high  SCL,  obtained  with  minimuin 
spatial  zonings  of  0.02  and  O.OO.S  m.  The  lield.-;  ;i.re 
evaluated  at  the  first  axial  grid  position  above  the 
emitting  surface.  See  Figure  I)-3  for  electric  field 
comparison  at  same  positions  for  tlie  two  cases. 


I'i^iuro  n :hI. 


Muj^iift  ic  fiold  time  liistorics  at  side  (if  eyliiuler 
. 4 m frcini  top  for  high  S(  1.,  obtained  with  minimum 
spatial  zoninj>s  of  0.02  and  O.OOS  m 


APPENDIX  E 

USER'S  MANUAL  FOR  ABORC  AND  PERIPHERAL  COMPUTER  CODES 


E-1.  INTRODUCTION 

This  appendix  contains  a user's  manual  for  ABORC  and  for  peripheral 
codes  PLOTALL  and  MOVIE.  Information  relating  to  specific  computer  require- 
ments for  the  system  is  given,  a flow  chart  is  presented,  inputs  are 
described,  and  a sample  calculation  is  illustrated.  The  particulars  of 
generating  overlaid  plots  and  creating  computer  MOVIES  are  also  spelled 
out  in  "cookbook"  fashion,  with  examples  provided.  The  PLOTALL  and  MOVIE 
codes  are  described  in  some  detail  because  they  can  be  used  to  operate  on 
data  from  other  electron  emission  and  SGEMP,  and  because  they  are  documented 
in  no  other  place. 


E-2.  ABORC  COMPUTER  REQUIREMENTS 

ABORC  is  a FORTRAN  IV  computer  program  operational  on  the  CDC  7600 
system.  The  code  consists  of  about  3500  cards,  including  all  FORTRAN 
programs  used.  No  machine  language  routines  are  employed.  The  code  cur- 
rently is  dimensioned  for  about  160,000  (octal)  words  of  fast  memory  stor- 
age and  about  720,000  (octal)  words  of  large-core  storage.  One  fast- 
access  file  is  required  during  execution,  although  two  are  used  if  the 
movie  option  is  invoked  by  the  programmer.  Run  time  of  the  cede  varies 
from  1 to  .30  CPU  minute.^  on  the  CDC  7600  computer,  depending  on  problem 
condi tions . 

An  automatic  warn  time  is  programmed  into  ABORC  so  that  sufficient 
time  is  allowed  to  save  plot  files  and  perform  graphics  functions  even 
if  the  computer  run  time  exceeds  the  maximum  amount  requested.  The  F'!N 
compiler  is  used  with  extended  and  preset-to-zero  core  options. 

L-.3.  ABORC  Fl.OW  CHART 

A brief  flow  chart  of  the  code  system  is  shown  in  Figure  L-1.  Note 
that  the  electron  emission  information  must  come  from  a source  external 
to  ABORC. 
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E-4.  ABORC  INPUT  DESCRIPTIONS 


The  problem  being  considered  is  defined  by  geometry,  emission  char- 
acteristics, and  conductivity  characteristics.  The  two-dimensional  geom- 
etry is  defined  by  specifying  annular  conductivity  zones  which  are 
described  by  minimum  and  maximum  axial  and  radial  coordinates.  The  total 
geometry  considered  by  the  problem  is  the  sum  of  all  defined  annular  zones. 
The  electrons  emitted  from  the  body  give  the  driving  function  for  the 
problem.  Emission  characteristics  include  axial  and  radial  bounds  of  each 
emission  zone,  energy  spectra,  angular  distributions,  and  emission 
intensities . 

The  method  of  describing  these  physical  properties  to  the  code  is 
outlined  below  in  the  input  card  description.  Variable  names  and  their 
meanings  are  given.  Data  formats  are  provided  if  required.  Name- list 
variables  are  i.iput  in  free  form  according  to  the  variable  type  (real  or 
integer),  with  commas  separating  them.  Column  1 must  remain  blank,  and 
tile  list  begins  with  b$DEFINE  and  ends  with  b$END  (where  b indicates  a 
blank  space) . 

The  word  DEBUG  is  used  at  the  beginning  of  the  description  of  those 
variables  whose  primary  function  is  use  in  debugging  problems.  Under  nor- 
mal circumstances,  the  default  values  of  these  inputs  arc  adequate  and 
the  variables  can  be  ignored.  The  word  liUIT  at  the  beginning  of  a vari- 
able description  mean.‘^  that  this  variable  is  used  in  editing  results  of 
the  calculation.s. 

The  names  in  parenthc.ses  are  the  variable  names  used  internally  in 
the  code  in  cases  where  the  internal  variable  name  differs  I'rom  the  injuit 
variable  name.  If  no  default  is  given  for  a variable,  that  variable  rnu^t 
be  input.  If  a v.iriable  has  a maximum  value  or  number  of  values,  the  des- 
ignation MAX  = or  MAX  NO.  + is  given.  A convenient  summary  of  the  minimum 
and  nr.’.ximum  values  is  found  in  Talile  li- 1 . 
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Table  E-l 

MINIMUM  AND  MAXIMUM  VALUES  OF  ABORC  VARIABLES 


Min  No. 

Min  No. 

Max  No. 

Input  Variables 


Zone  number  (each  direction) 

3 

100 

Time  steps 

- 

- 

Emission  zones 

1 

20 

Energy  distributions 

1 

20 

Time  histories 

1 

20 

Angular  distributions 

1 

20 

Conducting  annuli 

0 

100 

2-D  prints 

0 

30 

Mini-prints 

- 

- 

Plot  quantities 

1 

40 

CalCLilational  Variables 

Total  number  of  particles  emitted 
in  a given  time  step 

0 

1000 

Number  of  particles  being  followed 
in  a given  time  step 

0 

17,300 

Number  of  point  pairs  per 
distribut i on 

1 

50 

Total  allowed  values  in  all  point 
pair  tables  for  emission 

- 

^2000 

Total  allowed  point  pairs  in  time 
history  plots  (counting  all  curves 
which  are  overlaid) 

~> 

4000 

Km 


Certain  terms  peculiar  to  this  code  are  used  in  the  input  card 
descriptions  repeatedly.  Definitions  of  the  terms  are  listed  below. 


Mini-print 


2-D  print 


c 

Ax 


Card  1 C8A10) 
TITLE 

Card  2 14012) 
IOPT(1)=0 
= 1 


I0PT(2) 


I0PT(3)=+1 
= -l 

IOPT(4)=0 
= 1 


IOPT(4)=0 
= 1 


10PT(5) 


A short  printout  of  certain  problem  parameters  plus 
requested  plot  quantities  can  be  requested  at  times 
independent  of  large  "2-D  prints"  (see  below) . Auto- 
matically printed  at  2-D  print  times. 

Spatial  distribution  printouts  of  fields,  currents, 
and  charge  densities  (very  large). 

g 

Speed  of  light  (.3  x 10  m/sec). 

Used  to  designate  both  Ar  and  Az,  the  spatial  zone 
increments . 


The  driver  of  the  problem  is  the  particle  emitter. 

The  driver  is  an  analytic  current  specification  up 
the  axis  (see  subroutine  TESTJ) ; used  for  debug 
purposes . 

Default:  lOPT(l)  = 0 

DEBUG.  Particle  motion  printout.  Print  fir.st  and 
last  I0PT(2)  particles. 

Default:  I0PT(2)  = 0 

Plot  all  input  distributions. 

Suppress  input  distributions. 

Default:  I0PT(3)  = 1 

Random  emission. 

Analytical  emission  (see  subroutine  TESTJ) ; used  for 
debug  purposes. 

Default:  I0PT(4)  = 0 

Random  emission 

Analytical  emission  (see  subroutine  TESTJ);  used  for 
debug  purposes. 

Default:  10PT(4)  = 0 

DEBUG.  Print  first  and  last  lOPi (5)  emission  particles 
every  emission  step. 
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10PT(6)>0 

I0PT(7) 

IOPT(8j=0 
= 1 
= 2 

lOPT(lO) 

(NESCAP) 

lOPT(ll) 

IOPT(i;5) 

IOPT(i4)- 

10PT(16) 

IOPT(20) 

I0PT(21}=1 

IOPT(27) 

I OPT (37) 


DEBUG.  Statistics  check  plots  of  emitted  particle 
energy  spectra. 

Default:  IOPT(6)  = 0 (no  print) 

Number  of  CP  seconds  to  leave  at  end  of  run  for  plots, 
files,  etc. 

Default:  I0PT(7)  = 20 

Do  not  check  particle  emission  statistics. 

Check  particle  statistics. 

Check  and  print  debug  print  of  particle  statistics. 
Default:  IOPT(8)  = 0 

Number  of  chancer,  a particle  has  to  emerge  from  a non- 
zero conductivity  region. 

Default:  lOPT(lO)  = 5 

Beginning  random  number  selector.  Generator  is  called 
lOPT(ll)  times  before  start  of  a calculation. 

Default:  lOPT(ll)  = 0 

Input  Check.  Plots  body  shape  if  >-l.  Currently  not 
operationa.  Set  to  -1. 

Array  index  limits  used  with  I0PT(13)  for  limits  in 
axial  and  radial  directions.  Permits  expansion  of 
object . 

Default:  Problem  index  limits 

Debug.  Debug  print  in  function  converting  real  space 
to  generalized  coordinates  if  >0. 

Default;  IOPT(20)  = 0 
Debug.  Print  in  EMITTER. 

Default:  IOPT(21)  = 0 (no  print) 

Write  j)lot  file  information  eve  y iiil’r(27)  light  time 
steps. 

Default:  IOPT(27)  = 1 

MOVIE  file.  Write  particle  information  to  file  TAPE20 
every  TOPT(32j  {)articlc  time  steps  if  >0. 

Default:  IOPT(32)  = 0 
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I0PT(33;  EDIT.  Print  every  IOPT(3.3)  axial  zones  outside  region 

with  axial  zone  indices  I0PT(34)  to  IOPT(35)  in  2-D  prints. 
Default:  I0PT(33)  = 1 

I0PT(34)  EDIT.  Print  every  axial  zone  with  index  between  IOPT(34) 
and  10PT(35),  inclusive,  in  2-D  prints. 

Default:  I0PT(34)  = 1 

IOPT(35)  See  I0PT(33),  I0PTC34) 

Default:  Niimber  of  axial  zone  boundaries 
I0PT(36-38)  Like  I0PT(33-35)  only  for  radial  zones. 

I0PT(39)  EDIT.  Number  of  significant  figures  desired  in  2-D  prints 

for  all  quantities  except  electric  fields.  Also  gi/es  the 

range  of  the  variable  to  be  printed  out  “ i.e.,  from  peak 
- r IOPTf391 1 

value  down  to  10  times  peak  value.  Useful  in 

limiting  printout  size.  Use  -1  for  ElO  format. 

Range:  1 < I0PT(39)  < 7. 

Default:  10PT(39)  = 4 

All  of  the  following  variables  are  in  NAMELIST  DEFINE: 

DTN  Light  time  step  (sec).  Code  adjusts  to  meet  stability 

criterion  if  necessary. 

Default;  DTN  = Ax  . //2  c x 0.95,  but  then  DTN  is 
min 

decreased  to  the  nearest  integer  divisor  of  the  par- 
ticle time  step. 

DTPART  Particle  time  step  (sec). 

TMAX  Maximum  simulation  time  (sec). 

DTPRNT  Time  increment  for  mini-prints  (sec). 

Default  is  DTPR  values. 

DTPR  Times  at  which  to  j)rint  2-D  prints  (sec)  . 

MAX  NO  = 20. 

2-D  print  at  last  time  step  is  automatic. 

DTMOVl;  Time  increment  for  plotting  i)article  positions  on  the 

])rinrer-plotter  (sec).  If  zero,  no  plots. 

PlAJliNC  Multiplier  of  all  emission  current  levels  from  all  zones. 

Dimensionless . 

Default;  1 


i 
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PG  Axial  zone  boundaries  read  in  in  ascending  older  fm) . 

MIN  NO  = 3 
MAX  NO  = 100 

RG  Radial  zone  boundaries  read  in  in  ascending  order  (m) . 

MIN  NO  = 3 
MAX  NO  = 100 

The  variables  from  here  to  the  end  of  the  NAMELIST  DEFINE  are  read  in 
only  if  electron  backscattering  is  being  included  in  the  calculation. 

REFQ  Fraction  of  charge  on  particles  reflected  off  conductor 

walls  each  collision.  Dimensionless. 

Default:  0 

EFRAC  Kinetic  energy  per  unit  charge  contained  by  backscattered 

electi'ons  as  a fraction  of  incident  energy.  Dimensionless. 
Default:  0.9 

ZBOT  Bottom  position  of  cylindrical  region  from  which  to  back- 

★ 

scatter  electrons  (m) . 

RSIDE  Outer  radius  of  cylindrical  region  from  which  to  back- 

scatter  electrons  (m). 

QABTF  Fraction  defining  range  of  charge  particles  which  are 

continued  in  the  calculation.  All  particles  with  charge 
loss  than  QABTF  times  maximum  particle  charge  are 
e 1 i ininated . 

Default:  0 

$F,ND  NAMELIST  DEFINE 


★ 

Particles  will  backscatter  only  upon  striking  conduction  regions  in 
such  a way  that  no  grid  point  is  having  a contribution  to  its  current  den- 
sity due  to  that  particle  and  if  tiie  particle  is  entering  the  region  defined 
by  ZBOT,  ZTOP,  and  RSIDli. 
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Comment  Card 


Card  Type  1 
After  DEFINE 

Card  Type  2 
After  DEFINE 


Card  Type  3 
After  DEFINE 


Col. 

1-80 


1-10 

11-20 

21-30 

31-40 

41-50 

51-60 


61-80 


Specification  of  conduction  regions.  ElO  format. 

MAX  NO  = 100 

MIN  NO  = 0 fi.e.,  no  conduction  regions  are  necessary) 
SIGIM  Cower  axial  position  of  conducting  annulus  (m) 
SIGIP  Upper  axial  position  of  conducting  annulus  (m) 
SIG2M  Inner  radial  position  of  conducting  annulus 
SIG2P  Outer  radial  position  of  conducting  annulus. 

SIG  Conductivity  value  assigned  (mho/m) 

DEFAULT:  1 . 0E+9*eQ/At . 

FD  Read  in  non- zero  value  if  another  annulus  is 
to  follow;  otherwise,  blank. 

Conurient  field 


1-80  Comment  card 


Ni  I'E : The  rest  of  the  code  inputs  are  free-form.  Columns  1-80  may  be  used 
and  blanks  are  unnecessary.  Quantities  may  be  listed  in  any  order  witlvin 
the  division  labeled  EMISSION  CHARACTERISTICS. 


Card  Type  4 

After  DEFINE  IICZ^.R^),  ER(Z2,R2),  EZCZ^.R^),  ■ • • 

This  card  sets  quantities  which  arc  to  be  written  onto 
the  plot  file  for  use  by  the  plot  routines.  The  quanti- 
ties that  may  be  written  to  the  file  a:e: 

Plot  Type 


1 

H 

Magnet 

ic  field 

(amp/m) 

2 

EZ 

Axial 

electric 

field  ( 

volt/ m) 

3 

ER 

Radial 

elcctri 

c field 

(volt/m) 

4 

JZ 

Ax  in  1 

current 

density 

(amp/m" ) 

5 

JR 

Radial 

current 

density 

(amp/m^ 

) 

6 

KZ 

Axial 

surface 

current 

densi ty 

(arnp/m) 

7 

KR 

Rad i al 

surface 

cii  rrent 

densi  ty 

(amp/m) 

8 

lEMlT 

Total 

emi ssi on 

current 

(amp) 

9 

WEMIT 

Tot:'l 

cumulat i 

ve  em i s s 

ion  ener, 

gy  (joule) 

Plot  Type 
10 
11 
12 

13 

14 

15 


WKINETIC  Total  kinetic  energy  of  electrons  (joule) 

WFIELDS  Total  energy  in  fields  (joules) 

WLEAVE  Total  cumulative  kinetic  energy  leaving  system  (joule) 
\T4AX  Maximum  potential  anywhere  along  axis  (volt) 

VMIN  Minimum  potential  anywhere  along  axis  (volt) 

I LEAVE  Total  current  leaving  system  (amp) 

Plots  1-7  are  input  as  H(z,r),  etc. 


Plots  8-15  are  input  as  lEMlT,  etc.,  with  no  position 
specification. 

To  finish  the  list,  put  END  after  the  last  plot. 


• Commas  are  net  necessary  at  the  end  of  a card. 

• At  least  one  plot  quantity  must  be  specified. 

• As  many  cards  can  be  used  as  needed  up  to  the  specification 
of  40  plot  quantities. 

The  rest  of  the  inputs  specify  the  electron  emission  parameters  of  the 


problem. 

Card  Typo  5 

After  DEl-INE  liMISSION  CHARACTERIS  TICS 


Tills  card  simply  informs  tin  ode  that  emission 
information  follows.  Cards  within  this  grouji  are 
not  numbered  because  the  individual  distributions 
can  be  input  in  any  order  within  the  groiqi. 


h-4.1  Emission  Intensity  nj^l’oint  Pair  Table 

The  emission  intensity  format  is  used  to  specif)'  the  time  histor\'  for 
the  emission  pulses;  n is  the  intensity  numbei',  of  which  up  to  Jd  are  allow- 
able, and  POINT  I’AIR  I'ABLE  is  the  time  history  in  amp/iir  versus  time  in  sec. 
El  may  be  usetl  to  tibbreviate  emission  intensit)'. 

) 

Tor  examjile,  ;i  JO-nsec  FKIIM  triangular  pulse  with,  a peak  of  D.l  anip/iif 
would  be  injiut  : 


FMI‘^Sl(dJ  IIJILNS  ilY  ' 


i’Oi  . 


The  t;d)le  may  contain  up  to  50  point  pairs,  allowing  high  resolut  ion 
to  fairly  complex  emission  time  histories,  and  ma\'  Iv  coat  iiiued  on  as  iiians’ 
cards  as  desired. 


Id’F 


An  abbreviated  form  is  available  for  specifying  emission  pulse  shapes 
2 

of  the  form  f(t)  = sin  [ (7f/2)  (t/r)  ] , where  t is  the  rise  time  and  FWMM  of 
the  pulse.  Simply  specify  Eln  = SIN2  (T)  where  n is  the  distribution  num- 
ber set  by  the  user  and  T is  the  rise  time  in  seconds.  Note:  Abbrevia- 
tions listed  in  Table  E-2  may  be  used  here. 

Table  E-2 

FLOATING  POINT  NOJLTIPLIERS  FOR 
ABBREVIATED  ABORC  INPUTS 

Letter  Factor  Power  of  10 

P -12 

N -9 

U -6 

K +5 

M +6 

The  letter  factor  can  be  used  to  replace 
the  power  of  10  desired  on  ini>ut  cards. 

i:-4.2  Energy  Distribution  n (bins  = i)  = Point  Pair^  ljUile 

Similarly,  the  energy  distribution  card  allows  t\)i'  arbitrary  speci- 
fication of  the  (.'mission  energy  spectrum.  Again,  n is  the  distribution 
number,  and  uji  to  20  are  allowed.  (BINS  - i|  is  optional.  liD  may  be 
used  to  ab!)  rev  i :it  ('  energy  d i st  r i but  iini . i speciiies  tint  tlu'i'c  shall  be 
i jnirticU's  emitteii  ])er  particU'  t i mi'  stej)  pi'r  emission  I’.oiie  for  the  ilis- 
tribution.  If  the  (BINS  = i)  input  is  omitted,  ^ he  default  v:i  lue  of  i is 
■1  . I’he  POlN'f  PAIR  TABI.l.  is  as  belOri',  with  I'el.itivi'  fri'ipieiKV  |numbi-i/ 
unit  energy  versus  eiu'i'g,)’  I eV ) | as  the  parameters.  fdr  e.xampli',  the  spi'c- 
t rum  of  Figure  l.-l  migh.t  be  input: 

1 Nl  la  a 1 ' I ' . M-  I I ni:  ’ _ ‘a  , . -I  , i a _ . a _ ; i _ 1 _ 

' 'I  I ' e 1 ‘ . ’ I , ' ) 

Note  the  use  of  the  "k"  wliich  is  intirpreted  a:;  Id  . '.arioiis  ether  letters 
ma\'  hi'  used  with  float  lug  point  numbers  in  the  emission  sect  ion;  ! he\'  ari' 
eiuimerat  i'll  in  Tahli  12.  Ilii'  t'xampli'  speitrum  is  giieii  lei  riha  1 i ed  to  a 


5 


2) 


input  to  ABOR 


peak  of  1;  however,  the  emitter  normalizes  all  distribution  functions,  so 
the  emission  current  density  has  the  value  specified  by  the  emission  inten- 
sity versus  time  curve  multiplied  by  the  FLUENCE  and  TIMES  factors.  The 
FLUENCE  parameter  has  already  been  discussed;  the  TIMES  factor  is  discussed 
below.  Desired  relative  height  for  all  distributions  is  all  the  user  need 
consider. 

ANGLE  DISTRIBUTION  n = POINT  PAIR  TABLE 

The  angle  distribution  card  is  shown  above.  AD  may  be  used  to  abbre- 
viate angle  distribution.  The  parameters  for  the  point  pair  table  in  this 
instance  are  the  number  of  particles  per  unit  solid  angle  versus  angle  in 
radians.  For  example,  the  distribution  of  Figure  E-2  is  given  by: 

ANGLE  DISTRIBUTION  3 - 0,0.6,0.75,1.1,0.45,1.5708,0 

Polar  angles  are  employed  in  the  coding  for  emission  specification. 

The  angle  on  the  given  sample  is  the  angle  from  the  emitting  surface  nor- 
mal. The  azimuthal  angle  distribution  is  tyTically  specified  uniform 
from  0 to  2tt  for  two-dimensional  calculations.  Again,  up  to  .SO  point 
pairs  may  be  used. 

The  emission  energy  spectra,  angular  disti ibutions , and  pulse  shapes 
are  specified  by  the  above  cards  input  in  any  order.  The  last  section  of 
inputs  draws  all  this  information  together. 

E . 4 . 5 Emission  Zones 

The  emission  zones  card  simply  specifies  that  emission  zone  informa- 
tion follow.  The  emission  zones  are  indicated  via: 

(z^,r^)  TO  I NT  i DELAY  x TIMES  y,  ED  j,  AD  k,  ADQ3  £. 

^l’^l’^2’^2  coordinates  of  the  zone.  i is  the  emission  intensity 

time  history  number  reference  to  be  used  for  the  zone,  delayed  by  x seconds 
and  multiplied  by  y.  The  DEIAY  and  TIMES  factors  are  optional,  j is  the 
energy  distribution  number  to  be  used  for  the  zone,  k is  the  angle  distri- 
bution to  be  used  for  determining  the  direction  of  each  particle  relative 
to  the  surface  normal,  and  £ is  the  distribution  reference  for  the  emission 
elecTron  direction  azimuthal  angle  about  the  surface  normal.  Up  to  20  emis- 
sion zones  can  be  specified. 
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For  example,  suppose  we  desire  an  emission  zone  position  on  the  top 
of  a cylinder,  from  the  axis  of  0.2-m  radius  and  at  vertical  position 
(for  a particular  problem)  0.5  m.  We  shall  reference  the  previous  distri- 
butions and  input 

(.3,0)  TO  (.5,. 2)  INT  5 TlMhS  10,  £07,  AD3,  ADQ3  2 

where  we  have  multiplied  the  current  density  by  10.  We  may  also  want  to 
emit  from  the  side  cf  the  cylinder  (e.g.,  radius  =0.2  m)  but  with  reduced 
intensity  and  delayed  by  a couple  of  nanoseconds: 

TO  (.3,.  2)  INT  5 DELAY  2E;-9  TIMES  .1,  ED7 , AD3,  AD03  2 

where  and  r^  for  this  emission  zone  are  obtained  from  the  previous 
zone's  ^■2’^2' 

Card  16 
END  OF  INPUT 

terminates  the  input  processing. 

E-4  4 Final  Remarks 

All  cards  may  be  continued  on  the  next  card  except  for  zone  specifi- 
cations. Also,  blanks  make  no  difference  because  they  are  deleted  imme- 
diately by  the  coding  upon  reading  the  cards.  More  than  one  energy, 
angular,  or  time  distribution  can  be  emitted  fiom  a given  spatial  zone 
by  simply  respecifying  the  same  zone  coordinates  on  an  additional  emission 
zone  card. 

The  follov'ing  information  is  helpiul  in  using  the  code. 

• Angle  distributions:  polar  0 - it/2 

azimuthal  0 - Ttt 

• Plot  file:  tape  7 

• Movie  file:  tape  20 

• Kinetic  energy  print  is  every  particle  in  system  including 
first  time  step  particles 

• Zone  indices  increase  in  direction  of  increasing  z,  or  r 


e-5  sample  problem 


A sample  ARORC  calculation  is  illustrated  in  this  section.  Code  inputs 
and  selected  outputs  are  listed  to  summarize  previously  defined  variables  in 
a test  problem.  This  case  was  selected  as  a user's  illustration  for  code 
input  and  output  interpretation  and  to  provide  a benchmark  test  problem  for 
comparison  when  conversion  to  a new  computer  system  is  necessitated. 

The  sample  problem  has  the  following  characteristics. 


Geometry 
Outer  boundary 

Emission 


Pulse  shape 
Pulse  length 
Grid  characteristics 
Number  of  axial  zones 
Numt)er  of  radial  zones 
Number  of  emission  zones 
Number  of  energy  bins 
Particle  time  step 


Minimum  zone  size 
Clear  time 
Plots 


A cylinder  3 m high  x 3 m in  diameter 

Cylinder  is  in  an  outer  cylinder  30 
m high  by  33  m in  diameter 

146/m^  emitted  with  a cos  6 polar 
angular  distribution  and  isotropic 
azimuthal  angular  distribution  from 
the  top  and  half  the  side.  The 
electron  energy  spectrum  is  from 
the  QUICKE2  code. 

- 7 

sin"'^- 

17-nscc  rise  time  and  FWHM 


64 

58 

5 

3 

0.2  nsec 


Maximum  time  of  calculation  50  nsec 


0.02  m 
100  nsec 

At  various  points  on  the  top,  side, 
and  bottom  of  the  inner  cylinder,  both 
the  normal  electric  field  and  the  sur- 
face current  densities  are  monitored. 
Also,  oti'er  heliiful  quantities  are 
specified . 


E-5.1  1 


n i scuss i on 


The  input  cards  to  model  the  problem  mentioned  above  are  shown  in 
l-'igure  L-3.  The  plot  code  inputs  are  also  shown  after  tlie  end- of- ret  ord 
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T + .5S,DT.2,0J5,  DZ,02,5  £»»,JK5F;> 
OTPART*,?  f-q,DTpMNT.i .t-R, 


8HM»T,8S  ,H>'<4,TP1  7N,  TC:OON,3IO, 

• 1 

SDtUNt  so.t.-'), 

OfPR»20,E»'J,«0,t-‘),  60, ‘-P, 

ELUETiC*!  ('“■ta. 
lDPT(lJ)«i,^0,«3,•5,0.Sn, 

pjB  *1  s,  ,•!«,(•  s s, » “tg,  f 0 1 ••p, , ,s,  »e ,,  *7 ,8i  "7,, 

•6,Si“6,  ,*6#6,^6,2,"3,6#*3,6,*i,  ,*a,6i**? ,a,*l  ,6,  "1  tU#»l , #*,6, 

,32,.«,,S1  ,,63,,  76,.8<.,  1,1,  1,3,1. 5S,  1.0, 

a,OS,2,3,^.65,3.,3,!>,‘(,,si,!»,S,,^,5,6,,7,,8.,‘>,,10,,ll,,U,,13,,l<i,,l5,> 

«(;■  0.,. 2, .3-5, .A?, .55, .65, ,73, ,85, .95, 1.03,  1,15, 
l,25,l,33,l,u,l,«5,l.«8,l,5,l,52,1.50,1.56, 

1 , 59, 1,615, 1 ,698, 1,69, 1 ,7u5, 1,82, 1 ,89.2,, 15,?, 1, 

2, 5, i?,  7, 2, 9, 3,1  ,3. 3,3. 5, 3, 7, 3. 9, <1,1, 9.3, 

0.5, 9, 8, 5, 2, 5, 6, 6,, 6, 9, 6. 8, 7, 2, 7. 6, 8, 5,9,5, 10. 5,11.5,1?, 5, 13. 5, 19, 5,15, 5, 
16,5, 

SEND 

conductors 

«3,  1,5 

PLOTS  (07,005) 

KR(0.,,75),tZ(,005,,75),5Z(-,5,l,50  ),  ER{-  5,1,505) 

KZ(.l ,5, 1 ,50) , tR(-l,5,l,505) 

82t-2,5i  1 .50) » EH{-2,5.  1 .505’,hR(.3,?  .,75).  £Z(»3,e  - .75) 

VMAX,VMIN,If NtT.ILtAVE 

E2(,03,.75),£Z(,05,.75),tZ(.07,,75),EZ(,ll,.75),  1Z(.17,.75),  EZ(  ,25, , 75). 

EZ(  ,5.  .75)  .».KINETIC.*“HtLDS,f.NO 
EMISSION  CMAMACTf.HIST  ICS 

AD1»0,.0,,  .0675, ,1736,  ,1795, ,3020,  .2elB,.50,  .3991, ,6028,  .9363, ,766, 

0,5256, ,866,  ,6109, ,9397,  ,6981, ,9898,  ,7859,1,,  .0727, ,9898,  ,9599, ,9597, 

1 . 097, ,866,  1,139, ,766,  1,222, .6928,  1,509, .5,  1,396, ,592 
AD2aO, 1,6.28,1, 
tIl»8*N2(  17N) 

EDI  (8IN3*3  )»  l,00f*03,  l,32F*lb,  1.20t*05,  l,33t*lb,  l,60Et03,  l.oOEtlb, 

2,00f.05,  2,69f*15,  2,90E*05,  l,93t»l5,  2,80E+03,  1,251*15,  3,20E*03,  8,C8E*19, 

3,60t*03,  5.70E*19,  9,00t*03,  3,96t*19,  v:,50E*03,  2,10E*19,  5,OOE*03,  1,27E*19, 

5.50E  + 03,  7,66t.*i5,  6,00t*03,  9,63F*13,  6,50t«CJ,  2,91E*13,  7,0)E*03,  l,75t*l3, 

7,50£*03,  9. 961*12,  8,001^03,  5,70t*12,  9.00f*05,  1,981*12,  1,00E*09,  6,0et+ll, 

1,101*09,  1,091*11,  1.20E*09,  1,3*1*10 

EMISSION  ZUNIS 

iO.,0,)  TO  (0,,l,5)INTl  ,tDl ,A01,AD03  2 

TO(  -.38,1,5)  INTI  OILAT  ,6iN  TIMESI,  ,EU1,  ADI, 

TOC  -,75,1,5)  INTI  DELfV  1,9N  TIM131,  .EOl,  ADI, 

T0(-1,13,1,5)  inti  OILAT  3, IN  TIMisi,  ,£0), 

TOl  -1,5,  1,5)  INTI  niLAT  9.9N  TihOSl. 

END  OF  INPUT 
*7At 
BBMAT 
1 


ADQ3  2 
ADU3  2 
A51,  ADQ3  2 
IDl,  ADI,  AOU5 


1 , -09 


BEST  AVAIUBLE  COPY 


I'ij^ure  1:-?. 


ABORC  input  foi-  nadelinp, 


s;imi)le  problem 
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(fiOR)  card.  This  plot  input  causes  the  current  ABORC  file  to  be  plotted 
without  comparing;  it  against  any  other  file  from  a previous  calculation. 

1 i - S . 2 0_u  r pu t l)e  s cript  i. o n 

Some  outputs  given  by  ABORC  are  listed  in  this  section.  Figure  E-4 
shows  quantities  printed  by  the  code  which  are  very  helpful  as  input  checks. 
Such  items  as  total  r.uml)er  of  zones,  time  steps,  plot  file  sizes,  and 
approximate  calculation  time  give  a quick  indicator  if  something  is  grossly 
wrong.  Minimum  zones,  time  steps,  and  approximate  clear  times  are  also 
printed . 

The  mini-print  is  shown  in  Figure  F.-5  for  a time  of  20  nsec  into  the 
simulation.  Definitions  of  the  variable  names  titling  the  numbers  are 
listed  in  Table  F.-3  in  the  order  they  appear  on  the  printout,  reading  from 
left  to  right.  Some  of  these  variables  are  most  useful  in  evaluating  cal- 
culational  numerical  quality. 

The  numbers  following  the  labeled  variables  are  the  values  of  the 
requested  plot  quantities  at  this  time  step.  Numerical  hash  often  makes 
the  magnetic  field  and  current  density  values  too  noisy  to  be  useful  in 
this  un-averaged  form,  but  other  quantities  are  usually  well  behaved. 
Basically,  these  numbers  provide  a back-up  in  case  the  plot  file  gets 
lost  for  some  reason. 

Spatial  distributions  of  fields  and  currents  at  20  nsec  are  shown 
in  Figure  i;-b.  Only  those  portions  closest  to  the  axis  are  listed  due 
to  the  voluminous  output.  The  axial  (z)  and  radial  (r)  values  are  listed 
to  the  left  and  above  the  tal)les  of  values.  In  the  integer-printout 
tables  (11,  ,1Z,  .)R,  and  CtlAkCli  IN  ZONF.l,  only  the  index  of  the  radial  zone 
is  given  due  to  lack  of  room.  The  code  inputs  must  be  referred  to  in 
these  cases  for  tlie  radial  position  in  real  space.  Whether  to  use  ;i 
zt)ne  boundary  or  center  for  a given  quantity  can  be  determined  from  iable 
F.-4.  The  designations  "Q1  1)  1 RliC  I'lON  , ” "q2  1)  1 Rl.CTlON  , ” "Q3  DIRICTION" 
refer  to  axial,  railial,  and  azimuth.'il  directions,  respectively. 


i,OMf  INPUT  CHECKS 

NO.  OF  PARTICLE  TIHE  STEPS  2Sn 

NO.  OK  FIELD  TIME  STEPS  12b0 

NO.  OF  ZONES  I6‘t8 

NO.  OF  TIME  STEPS  TO  BE  PUT  ON  PLOT  FILE  F25 

NO.  OF  WORDS  TO  BE  WRITTEN  TO  PLOT  FILE  1.4J7SE*0R 
NO.  OF  FIELD  CALCULATIONAL  STEPS<70NE  STEPS)  R.5£C0E»06 
TIME  TO  CALCULATE  FIELDS(7feOO  MIN)  C.BAOCt-Ol 
APPROX  COSTIDOLLARS)  1.B240Ea01 

APPROX  COST  PER  non  AVERAGE  PARTICLES  THIS  RUN  ( DOLLA  BS  > .T.;E0rt*01 
flEAR  times  are  less  than  2.001‘iE-07  1.10PH)-07  in  the  axial  and  RADI 


21  72  R2 

-1.5000E+01  l.SOnOftOl  1.6500E*U1 

THE  SMALLEST  ZONE  IN  THE  RADIAL  DIRECTION  IS  2.00000nE-02 
THE  SMALLEST  ZONE  IN  THE  AXIAL  DIRECTION  IS  ?,OOOOOOE-n2 

particle  time  S1EP=  2.00000E-10  SECONDS 

LIGHT  TIME  STEPr  R.OOnoOE-ll  SECONDS 

UPCATE  PARTICLE  POSITIONS  EVERY  S LIGHT  TIME  STEPS 

MAX  ALLOWABLE  TIME  STEP  FOR  STABILITT  IS  4.’)H144E-U  SECONDS 


F'igure  E-4. 


ABORC  printout  of  cpiantities  which  test  tht; 
inputs  for  ri asonableness 
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Figure  E-5.  Example  of  ABORC  mini-print 


Table  E-3 

MINI-PRINl’  VARIABLE  GLOSSARY 


NSTEP 

TIME 

NEMIT 

NPART 

I EM  IT 

QEMIT 

QINSYD 

QLEAV 

QCHEK 

ENERGY  EMIT 
ENERGY  KINETIC 
ENERGY  ABSORBED 
ENERGY  EFIELDS 

ENERGY  HE I ELDS 

ENERGY  TOTAL  FIELDS 

VMAX 

POS 


ECONSERVE 


Light  time  step  number 
Simulation  time  (sec) 

Number  of  particles  emitted  during  present  time  step 

Number  of  particles  being  followed  at  this  time 

Total  emission  current  at  this  time  (amp) 

Cumulative  charge  emitted  up  to  this  time  (coulomb) 

Total  particle  charge  representing  emitted  electrons 
at  this  time  (coulomb) 

Total  emission  electron  charge  striking  and  sticking 
to  material  surfaces  during  this  time  step  (coulomb) 

Charge  conservation  check;  not  entirely  correct  at 
present,  but  small  values  indicates  good  conserva- 
tion of  emission  charge 

Total  cumulative  kinetic  energy  emitted  up  to  this 
time  (joule) 

Total  kinetic  energy  of  all  emission  electrons  at 
this  time  (Joule) 

Total  cumulative  emission  electron  energy  absorbed 
by  materials  up  to  this  time  (joule) 

Total  energy  stored  in  electric  field  at  this  time 
(joule) 

Total  energy  stored  in  magnetic  field  at  this  time 
(joule) 

Total  energy  stored  in  E and  M fields  at  this  time 
(joule) 

Maximum  electric  potential  anywhere  along  the  axis 
relative  to  the  bottom  of  the  outer  can  at  this 
time  (volt) 

Position  of  maximum  potential  along  the  axis  rela- 
tive to  the  bottom  of  the  outer  can  at  this  time  (m) 

Energy  conservation  parameter:  fraction  of  energy 
emitted  minus  energy  in  fields  minus  kinetic  energy 
minus  energy  absorbed  to  total  energy  emitted  (dimen- 
sionless); small  value  ('vO)  indicates  conservation 

Energy  in  emission  electrons  and  fields  i:omparcd  to 
energy  emitted  minus  energy  absorbed  (dimensionless); 
small  value  ('^0)  indicates  conservation 


ECONSERVE  IN 


*WIT '^,’!'T>  ' ? V ■?  W/'^- -M  1 1»  p 1 rjT- n njr 


NLEAV 

I.EAVT 

NEMITT 

NSTEPP 

QLEAVT 

NABORT 

N ABORT! 

NREFT 

QREFl' 

IQABTT 

QABTT 

QMAX 

QMAXT 

QMIN 

QMINT 

V 


Table  E-3  (cont.) 


Number  of  particles  oeing  absorbed  by  material  during 
this  time  step 

Cumulative  number  of  particles  absorbed  up  to  this 
time 

Total  number  of  particles  emitted  up  to  this  time 
Particle  time  step  number 

Cumulative  emission  electron  charge  absorbed  by  con- 
duction surfaces  up  to  this  time  (coulomb) 

Number  of  particles  prematurely  aborted  during  this 
time  step;  see  lOPT(lO) 

Number  of  particle  re-emissions  of  electrons  occur- 
ring during  this  time  step  if  electron  backscattering 
is  being  treated 

Cumulative  number  of  re-emissions  of  electrons  occur- 
ring during  this  time  step  if  electron  backscattering 
is  being  treated 

Cumulative  amount  of  charge  re-emitted  up  to  this 
time  due  to  electron  backscattering  (coulomb) 

Cumulative  number  of  particles  aborted  up  to  this 
time  due  to  charge  becoming  less  than  QABTF*QMAX, 
where  QMAX  is  defined  below  and  QABTF  is  defined  in 
the  input  description  section 

Cumulative  amount  of  charge  aborted  up  to  this  time 
(coulomb) ; see  lOABTT  above 

Charge  on  the  smal Icst-magni tude  particle  being  fol- 
lowed during  this  time  step  (coulomb) 

Same  as  QMAX 

Charge  on  the  largest -magnitude  particle  being  fol- 
lowed during  this  time  step  (coulomb) 

Same  as  QMIN 

Maximum  velocity  of  all  particles  being  followed 
during  this  time  step  (m/sec) 
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Figure  E-6.  Spatial  distributions  of  fields,  currents,  and  charge  for  ABORC  calculations 
at  20  nsec.  This  printout  is  also  designated  "2-D  print"  elsewhere  in  this 
appendix. 
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Figure  E-6  fcont . 
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Figure  E-6  (cont. 
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BEST  AVAIUBLE  COPY 
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Table  E-4 

RELATIVE  POSITIONS  ON  THE  GRID  WHERE  FIELDS, 
CURRENTS,  AND  CHARGE  ARE  CALCULATED 


Quantity 


Axial  Radial 

Position^  Position^ 
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CHARGE 
IN  ZONE 
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b 
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c 
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b 

c 
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E-5 . 3 PLOTALL  Output  Description 

The  plot  code  output  is  shown  in  Figure  E-7  as  it  typically  appears  at 
the  end  of  an  ABORC  calculation.  This  example  is  for  the  case  of  plotting 
only  one  file  of  data  generated  by  a single  ABORC  run.  Discussion  on  over- 
laying additional  files  of  data  is  given  in  the  following  section. 

Notice  that  the  plot  code  inputs  are  printed,  and  then  most  of  the 
ABORC  inputs  are  given.  Those  quantities  reside  on  the  plot  files, 
enabling  the  user  to  check  an  old  file  to  determine  exactly  what  ABORC 
inputs  were  used  to  generate  it.  Also  notice  that  the  exact  plot  posi- 
tions are  listed  along  with  information  pertinent  to  file  and  ])lot 
man  ipulati on. 

Finally,  sample  plots  Ui'e  shown  in  Figure  E-8  for  tiie  normal  electric 
field  at  the  top  of  the  cylinder  and  the  surface  current  density  at  the 
middle  of  the  side.  The  curves  are  smoothed  over  1 nsec.  liu'  plot  tyjie 
and  exact  leiation  on  the  ABORC  grid  are  given  above  tlie  plot,  along  with 
;i  convenient  plot  number  and  the  ABORC  title  card. 
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E-5  GENERATION  OF  OVERLAID  PLOTS  ~ PLOTALL  CODE 

A FORTRAN  IV  computer  code  called  PLOTALL  is  part  |of  the  ABORC  com- 
puter system.  This  code  operates  on  ABORC  and  QI.IICKE2  (Ref.  1)  data  files 
to  produce  plots  of  specified  quantities.  As  many  as  three  such  files  can 
be  overlaid  on  the  same  grid  and  the  plots  directed  to  the  printer  or  to  a 
pen-plotter.  Computer  requirements  for  these  plots  are  minimal  for  printer 
plots  and,  of  course,  require  a pen-plotter  for  the  higher- resolution  plots. 
Also,  piors  can  be  made  on  CRT  or  ether  plotters  with  appropriate  software 
used  with  PLOTALL.  This  section  discusses  how  to  generate  overlaid  plots 
with  the  code.  Inputs  are  described  and  a sample  case  is  discussed.  Empha- 
sis is  given  here  to  an  example  of  overlaid  plot  files.  An  example  of  plot- 
ting a single  file  is  given  in  the  previous  section  in  conjunction  with  the 
ABORC  sample  problem. 

While  the  examples  given  in  this  appendix  are  for  ABORC  data  file 
manipulation,  the  capabilities  discussed  and  il lastrated  also  apply  to 
QUICKE2  data  files.  In  the  latter  case,  spectra  are  plotted  rather  than 
time  histories,  and  tyincally  no  data-smoothing  is  performed.  PLOTALL 
capabilities  such  as  comparing  curve  1 on  file  TAPE?  with  curve  2 on  file 
TAPES  are  still  the  same.  Most  plot  runs  are  dene  using  the  negative 
value  fox  the  "number  of  fixes"  specification  [lOPT(l)],  liowever,  which 
allows  inputting  only  tw  iata  earns  to  define  the  plots  once  the  proper 
data  files  are  loaded  on  the  system. 

An  important  function  performed  by  the  PLOTALL  code  is  that  of  curve- 
smoothing. This  optional  capability  is  often  used  in  attempts  to  filter 
out  erroneous  high-frequency  noise  in  time  history  data  from  ABORC  due  to 
particle  motion  through  finite  spatial  grids.  This  noise  is  a byproduct 
of  following  too  few  particles  to  permit  cancellation  of  statistical  fluc- 
tuations. The  smoothing  formula  employed  by  PLOTALL  is  equivalent  to  the 
t ime -averaging  of  a voltage  pulse  entering  a circuit  with  parallel  resis- 
tance and  capacitive  elements.  The  capacitor  voltage  responds  to  the 
incident  voltage  with  a characteristic  RC  time-  constant.  Similarly,  all 

^S.  H.  Rogers  and  A.  J.  Woods,  "QUICKE2'  .An  Analytical  Electron  Emis- 
sicn  Code,"  INTEL-RT  8141-026,  .June  15,  1976. 
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curves  treated  by  PLOTALL  can  be  smoothed  as  if  they  are  voltage  pulses, 
causing  a delayed  response.  The  characteristic  time  is  designated  as  the 
averaging  time  and  is  specified  to  the  code.  The  equation  for  the  aver- 
aging is 


df(t)  ^ 
dt 


f(t)  _ 


g(t) 


where 

g(t)  = input  function  to  be  time-averaged  smoothed, 

T = characteristic  smoothing  time, 
f = resulting  time- smoothed  function. 

E - 6 . 1 PLOTALL  Inputs 

Inputs  to  the  PLOTALL  code  are  described  here,  along  with  a few  pointers 
on  required  specifications  for  desired  outputs. 

In  the  input  card  descriptions  which  follow,  the  abbreviations  D and 
MAX  stand  for  default  and  maximum  allowed  value.  Notice  the  special 
instructions  at  the  far  right  of  the  figures  for  some  of  the  cards.  Addi- 
tional instructions  are  given  below,  along  with  Table  E-5,  which  lists 
some  constraints  on  input  and  calculation  variables. 


Note  on  De faults 

Card  types  4 and  5 must  be  input  in  groups  of  NPILE  cards  of  each 
type,  where  NPILE  can  be  from  1 to  5.  Values  on  cards  after  the  first 
one  in  each  set  are  assumed  to  be  the  value  on  the  first  card  unless  a 
value  appears  on  the  subsequent  card.  Therefore,  defaults  listed  here 
pertain  only  to  the  first  card  of  the  set. 


Additional  Notes  and  Definitions 


• PLOTALL  looks  for  data  on  files  TAPE.'/,  TAPES,  and  TAPE  . depend- 
ing on  number  of  files  being  overlaid. 

• Curves  on  the  same  file  can  be  overlayed  by  copying  the  data  to 
two  of  the  files  used  by  PLOTALI.  and  then  specifying  appropriate 
plot  pcs  it  ions  with  the  IPl’OS  cards  (card  4). 

• "I  i le"  is  defined  as  a block  of  curves  generated  by  one  ABORC  or 
one  QU1CKE2  calculation. 

” "['lot”  is  defined  as  the  figures  generated  in  which  one  to  three 
curves  of  the  separate  "files"  are  overlaid. 
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E-6.2  PLOTALL  Sample  Problem 

A sample  problem  is  considered  here.  The  example  is  to  overlay  three 
separate  ABORC  results  on  the  same  grid.  Each  file  is  copied  to  one  of  the 
three  used  by  PLOTALL,  TAPE7-9.  The  inputs  to  the  code  are  shown  in  Figure 
E-9.  The  spacings  indicate  that  blank  cards  have  been  used  for  specifica- 
tions for  data  on  TAPES  and  TAPE9.  The  defaults  of  plotting  every  curve  on 
every  file  and  averaging  over  the  same  time  constant  have  been  employed. 

An  example  output  of  the  computer  run  is  given  in  Figure  E-10.  Notice  the 
plot  title  at  the  top,  the  plot  number  given  by  the  code  for  convenience, 
the  averaging  time,  the  plot  type,  exact  position,  and  ABORC  run  titles 
defining  the  curves.  All  of  these  features  are  printed  automatically, 
maki  g identification  of  the  plots  very  positive.  Additional  features  of 
the  PLOTALL  code  are  described  in  the  ABORC  sample  problem  section. 


E-7  GENERATION  OF  ELECTRON  MOTION  MOVIES  - MOVIE  CODE 
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A FORTRAN  IV  computer  code  called  MOVIE  is  part  of  the  ABORC  system. 
This  code  operates  on  ABORC  and  DYNASPHERE  (Ref.  2)  data  files  to  produce 
movies  showing  the  motion  of  particles  representing  photo-electrons  emit- 
ted from  the  objcct(s).  The  code  requires  a file  of  data  (cither  a per- 
manent file  or  a tape)  produced  by  invoking  the  appropriate  option  in  the 
SGFMP  calculation  and  taking  the  simple  control  card  procedures  necessary 
for  saving  the  file.  This  section  details  to  theprogrammer  how  to  make  a 
movie  of  ABORC  jihoto-el  ect  rons  once  tlic  file  of  information  is  available. 
Generation  of  the  I'ile  is  cxiUained  above  in  the  ABORC  input  description 
section  or  in  the  DYNASPllIiRFi  user's  manual.  A brief  description  is  given 
below  of  steps  necessary  to  use  the  MOVIFi  code,  description  of  code  inputs, 
sample  inputs,  and  discussion  of  a sample  movie  frame. 

H-7.1  MOVIE.  Injjut  Rc‘.|i.i i reinent s 

ABOlU;  provides  all  the  particle  [losition  information  at  each  time 
step,  but  the  MOVlli  code  is  not  presently  set  up  to  automatically  JescriFie 
the  bod)'  outlines  from  the  ABORC  inputs.  Therefore,  certain  inputs  must 

"A.  ,1.  Woods,  "User's  Manual  for  the  DYNA.SPIIT.RI;  SCdiMl’  t.om[)Uter  Code," 
INTLL-RT  8141-029,  April  1970. 


LOT ALL  INPUTS 


ime  ^cr  each  data  point,  then  TSTART  and  TFIN  are  checked. 


Table  E-5 

MINIMUM  AND  MAXIMUM  ALLOWED  VALUES  OF  PLOTALL 
INPUT  AND  CALCULATIONAL  VARIABLES 


Quantity 

Minimum 

Value 

Allowed 

Maximum 
Value 
A1 lowed 

Number  of  different 
files  overlaid 

1 

3 

Number  of  different 
plots 

1 

40 

Number  of  point  pairs 
on  one  plot 

2 

4000 

hbmat,hlkj  sens 

3 


t ijiuri,'  li-D.  Samjilt'  PLUrALI-  input  card.  Ihose  cards 
indicate  that  three  files  of  data  are  to 
be  overlaid  and  avcraj^jed  over  2 nsec, 
riiere  are  three  blank  cards  after  the  "3" 
ami  two  blank  cards  after  tlie  "2.01)"  card. 


coooo 

tooc  DOOOOCCU 


Figure  E-10=  Sample  output  of  computer 


be  specified  to  MOVIE  if  an  outline  of  the  object (s)  is  desired  on  the 
frame.  Other  options  such  as  spatial/temporal  plots,  which  can  appear  on 
the  frame  along  with  the  particle  motion,  must  also  be  specified  to  the 
code.  These  capabilities  can  be  invoked  using  straightforward  procedures 
described  in  Table  E-6.  Obviously,  the  additional  steps  of  special  con- 
trol card  procedures  are  needed  to  direct  the  calculations  to  the  proper 
disposal  station  of  the  computer,  but  these  commands  are  entirely 
installation-specific  and  so  are  not  discussed  here.  The  capabilities  of 
specifying  time  histories  or  spatial  distributions  of  plot  quantities  along 
with  particle  trajectories  exist  only  on  ABORC  MOVIE  files.  DYNASPHERE 
movies  contain  only  particle  trajectors,  and  the  spherical  surfaces  are 
obtained  automatically  from  the  data  file.  All  of  the  inputs  described 
in  Table  E-6  are  still  defined,  but  some  are  only  dummies. 

E - 7 . 2 Sample  MOVIE  Code  Inputs 

A sample  of  MOVIE  code  inputs  is  shown  in  Figure  E-11.  Note  that 
this  input  does  not  produce  the  sample  frame  shown  in  Figure  E-12. 

E ••  7 . 3 Sample  MOVIE  Frame 

A sample  frame  generated  by  the  MOVIE  code  is  shown  in  Figure  E-12. 

The  axis  of  the  object  is  at  the  left  side  of  the  frame.  At  this  time  of 
28  nsec  (shown  at  bottom  center),  the  electrons  have  moved  far  away  from 
the  object.  This  case  was  for  low  fluence,  .-ui  no  space-charge- 1 imi  t ing 
occurred.  The  graph  on  the  right  shows  the  magnetic  I'ield  near  tlie  boom 
of  the  object.  Any  quantity  availal)l('  on  the  data  file  could  have  been 
specified  here.  ihe  large  title  at  the  top  would  appear  in  the  position 
of  the  first  title  card  shown  in  Figure  li- 1 1 , while  the  smaller  title 
under  the  graph  would  be  iiqmt  in  the  second  title  card  iiqniT  given  in 
the  example.  The  small  title  iiiidenie.ith  the  main  title  on  the  frame  is 
the  t i t U'  caixl  used  in  the  ABORT  run.  which  is  automat  icall>’  obtained 
from  the  MOVIE  I'lle  by  tlie  coding. 


Table  E--6 

INPUTS  FOR  MOVIE  CODE 


Card (s) 

Columns 

F’ormat 

Variable 

Description 

1 

1-32 

8A4 

f40VTIT(I) 

Movie  title  (placed  at  top  of  each 
frame) . 

2 

1-10 

110 

I PART 

>0  to  plot  trajectory  of  particles. 

50  no  particles  plotted. 

2 

11-20 

no 

NGRAPH 

Number  of  curves  to  be  plotted  in 
graph  if  IPOS  < 0;  number  of  points 
on  graph  if  IPOS  > 0;  NGRAPH  - 15. 
ABORC  only. 

2 

21-30 

no 

IPOS 

>0  plot  graph  quantities  versus 
position  (or  other  measure). 

50  plot  graph  quantities  versus 
time.  ABORC  only. 

3 

1-24 

6A4 

GRATIT(I) 

Graph  title  (placed  below  graph). 

ABORC  only. 

4 

1-80 

8110 

I GRAPH (I) 

File  numbers  of  data  to  be  plotted 
on  graph  (files  are  read  from  ABORC 
plot  tape);  1 to  NGRAPH  values  up 
to  3 allowed.  ABORC  only. 

(UhO  card 

S on  1 y 

if  Il’OS 

0) 

5 

1-80 

8E10. 3 

POS  ( 1 ) 

I’osition  (or  other  measure)  of 
graph  quantity  1 ; 1 to  NGRAPH  values. 

6 

1-10 

E 10.3 

TAV 

Time  to  average  quantities  over  (nsec) 
ABORC  only. 

() 

11-20 

lilO.  3 

I'GliOM 

Factor  to  reduce  bound  of  geometry  liy . 

7 

1-10 

1 10 

lRi;i> 

Number  of  times  to  repeal  each  movie 
frame. 

8 

1-10 

1 10 

NCi.osi: 

■0  if  close-up  of  part  of  the  body 
lies  i red  (give  only  the  geometry 
of  that  |)art) 

: 0 whole  body. 

8 

11-20 

1 10 

N.XY 

Total  number  of  points  to  define 
geomet  1')'  ol  hodv  (only  i n t e isec  t i on 
points  need  be  given);  ABORC  only. 

8 

2 1 30 

1 1 0 

NG 

Nnmbei'  v)f  groiq's  in  which  to  construct 
body  (points  in  each  group  are  con 
necteil);  1 Nii  ' 1;  .ABORC  oiil)'. 

s 

3 1 - 70 

!■  1 1 0 

N’I’l  1 ) 

Numher  ol  points  in  each  groi  a;  1 to 
.NG  values;  .ABORC  onl\'. 

') 

1 -80 

81  10.3 

7,r(  1 ) 

7 cool'll  in.it  i‘  of  jioints  det'ining  lioih' ; 

1 to  .N.XY  V, lilies;  ABORC  only. 

1 (1 

1-80 

81  1 ('.3 

P.IM  1 1 

K coorilinate  ol  points  definiiu;  hodv; 

1 to  :\XV  valuos;  ABOUC  oiiIn', 
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Figure  F-12.  Sample  movie  frame  produced  by  MCVll  (.ode 
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